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FOREWORD 


This report presents the second volume of a two-volume final 
report on a one-year study entitled "Analysis of Apollo Space- 
craft Parachutes." (The companion volume Is listed as Reference 
1. ) This study was performed by Northrop Ventura for NASA/MSC 
under Contract NAS-9-8131. Messrs. M. A. Sllvelra, K. Hinson 
and C. Eldred of NASA/MSC monitored and reviewed the study. 

This study, designated as Project 0111 at Northrop Ventura, 
was carried out with direction from the Systems Engineering 
Group under Mr. R. G. Iemm. Program direction was provided 
by Mr. T. W. Knacke of the Advanced Design Group, and the 
Project Engineer was Mr. F. E. Mickey of the Aerospace Landing 
Systems Project Office. 

Principle authors of the sections in the body of the report 
are as follows: 

Mr. W. M. Mullins, Section 3; 

Mr. D. T. Reynolds, Section 4; 

Dr. K. G. Lindh, Section 5; and 
Mr. M. R. Bottorff, Section 6. 

The authors greatfully acknowledge the contribution of Mr. 

R. L. Ranes in the development of the structural analysis 
methods presented in Section 3 and also the valuable assistance 
of Mr. A. J. Me Ewan, who prepared Section 5.5- 
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ABSTRACT 


Hie major factors in structural analysis of parachutes are re- 
viewed, including ultimate strength, design factors, safety 
factors, external loading, and internal loading of the individual 
structural members. A method for computing the internal load 
distribution for a parachute with a given applied riser load 
and canopy differential pressure distribution is presented, 
and comparisons of the predicted canopy shapes and Internal 
loadings with Apollo development test results are given, lhe 
analysis discussed above is modified to give better representation 
of ribbon parachutes with vertical ribbons, and a comparison of 
results of the two analyses is given. A study of data obtained 
in dynamic loading tests of pilot chute risers adds to the under- 
standing of the effects of high load onset rates which occur 
at main parachute canopy stretch. A study of methods for measuring 
parachute aerodynamic and internal loading is performed, and a 
plan is developed for a three-phase test program to acquire needed 
data. 
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SECTION 1.0 
INTRODUCTION 


Hie Apollo parachute landing system was designed, developed and 
qualified by Northrop Ventura during the period 1962-1968. In 
the normal course of this development, many flight tests were 
made, and extensive data on the performance of the Apollo space- 
craft parachutes were collected. Hiese data were used as the 
basis for developing the methods that have been used during the 
course of the flight test program for estimating loads and in 
maklr£ structural analyses for the three Apollo parachute as- 
semblies: the drogue* pilot and main parachutes. 

It was recognized that there would be substantial value in an 
analysis effort that would review all the flight tost data at 
one time. In particular, it was seen that an analysis effort 
at this time would be free of the day-to-day pressures associated 
with a development program, and that consequently It could up- 
grade the loads and stress analysis methods used for the Apollo 
parachutes in ways that could not have been done previously. Hie 
present study was therefore authorized with the objective of up- 
grading and Improving the loads, stress and performance pre- 
diction methods for the Apollo spacecraft parachutes. Also, in- 
cluded in this study are the three additional tasks: (1) developing 

ideas for a new theoretical approach to the parachute opening 
process, (2) developing ideas for new experimental-analytical 
techniques to improve the measurement of pressures, stresses and 
strains in inflight parachutes, and ( 3 ) conducting a computer 
study to explain the data obtained in certain riser dynamics tests 
performed at Northrop Ventura in 1967 . 
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The results of the study are published in two volumes as follows: 


and 


INVESTIGATION OF PREDICTION METHODS FOR THE LOADS 
AND STRESSES OF APOLLO TYPE SPACECRAFT PARACHUTES 

VOLUME I - LOADS 


INVESTIGATION OF PREDICTION METHODS FOR THE LOADS 
AND STRESSES OF APOLLO TYPE SPACECRAFT PARACHUTES 

VOLUME II - STRESSES 

The present volume Is VOLUME II - STRESSES. The companion 
volume is listed as Reference 1. 

1.1 INTRODUCTION TO VOLUME II 

Increases in the Apollo spacecraft weight without corresponding 
increases in recovery system weight and volume brought about a 
need for increased accuracy in the structural analysis of para- 
chutes. Improvements in the prediction of material and joint 
strengths were made by laboratory testing and by use of prob- 
ability theory. Improvements in the predictions of internal 
load distribution were accomplished by development of improved 
analytical methods. 

The internal load analysis method that was evolved during the 
Apollo development programs was a significant advance in the state 
of the art. However, the full potential of the method was not 
utilized because of the press of schedule. A large number of 
hand calculations were required to obtain a solution so that 
only a limited number of conditions could be investigated, even 
with a simplified structural model. Correlation of the analytical 
results with test data was also limited. 

The purpose of this study is to further improve the structural 
analysis methods. Use is made of the Apollo flight test data to 
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corroborate and refine the structural analysis methods. A digital 
computer program developed during this study is used to investi- 
gate the internal load distribution during the deployment process 
to a degree not previously possible. This analysis is x'urther 
extended to give more realistic modeling of ribbon parachutes. 

An investigation of pilot chute riser dynamics utilizing Apollo 
laboratory test data is also included. The need for additional 
information on the aerodynamic forces acting on parachutes 
during deployment is demonstrated, and a plan for acquiring this 
data is developed. 

Section 1.2 is a review of the literature on parachute structural 
analysis, starting with the work of Jones and Taylor published in 
1923. In Section 2 the major terms used in the structural analysis 
of parachutes are defined and discussed. 

The methods for determining internal loads used in this study are 
derived in Section 3, and flow diagrams of two computer programs, 
CANO and CANO 1 which implement the analyses are given. 

In Section 4 the analyses discussed above are applied to the 
Apollo drogue, pilot and main parachutes, and the results are 
compared with test data accumulated during the Apollo development 
and qualification test programs. Included in Section 4 are the 
following analyses: Drogue chute reefed — interal load distri- 

bution and comparison of shape with test photographs; Drogue chute 
disreefed — modeling sensitivity study, pressure distribution 
sensitivity study, canopy growth study, film analysis, failure 
analysis, and vertical riobon load analysis; Pilot parachute — 
interal loads analysis and film analysis; Main Parachute internal 
load analysis, failure analysis, pressure distribution and film 
analysis, stress-time history, and optimum weight calculations. 

In Section 5# a study of the dynamic loading of pilot parachute 
risers is given. Section 6 is a study of the techniques for 
measuring the parachute canopy differential pressure and internal 
loads. 
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1.2 LITERATURE SURVEY - STRESS ANALYSIS METHODS 

The stresses In a parac.iute canopy are dependent on the Inflated 
shape, which in turn depends on the stress distribution. This 
interaction requires that the solutions of canopy stresses and 
shape be obtained simultaneously. It is understandable, there- 
fore, that much of the effort in parachute stress analysis is 
concerned with canopy shapes. 

The earliest analytical studies of parachute canopy shapes were 

performed by the British investigators Taylor, Southwell, Griffiths 

Jones, and Williams at the Royal Aircraft Establishment in 1919. 

2 

This work was summarized by Jones in 1923 • The theoretical 
canopy shape (known as the Taylor Shape) that was derived has 
served as the basis for most theoretical work until recent times. 

The "Taylor Shape" is the shape that would be assumed by an axi- 
symmetric parachute with an infinite number of suspension lines 
and zero hoop stress. Taylor reasoned that this shape, which is a 
surface of revolution in the limit, would be the flattest possible 
and therefore the optimum profile. He then concluded that the 
excess material in the pleats or bulges could be removed without 
affecting the profile, thereby giving a parachute of minimum weight 
Although his reasoning at this point was erroneous, (the effect 
of the excess material in the bulges does not disappear as the 
number of suspension lines approaches infinity), the shape derived 
gives a good approximation of the radial tape profiles of con- 
ventional parachutes. 

Taylor also formulated the equations for equilibrium and continuity 
in a gore bulge and showed that his equation in general form is 
equivalent to the general membrane equation. 
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If the differential pressure is uniform over a canopy, the Taylor 
shape is defined by: 

2 

^ = sin 0 

a d 

where : 

r = the local half -diameter 

a = the maximum half -diameter 

= angle between the central axis 
and a normal to the canopy 
surface 

Tables of rectangular coodinates for the profile curve are given 
in Reference 2 for the constant pressure case and two nonuniform 
distributions expressed in functional form. 

Equations are also developed by the other authors for shapes in 
which the hoop stress is not zero, and tables of rectangular 
coordinates are given for cases in which the hoop tension is 
constant or varies as a function of r. Weight estimates are 
presented for parachutes constructed to these theoretical shapes 
and cloth stresses are given. Section 3 of Reference 2 gives 
results of wind tunnel tests of these parachute shapes, including 
pressure and force measurements. 

In 19 ; +2 Stevens and Johns^ considered the case of a parachute with 
a finite number of suspension lines which extended over the canopy. 
The cords were made shorter than the gore panels so that the cloth 
carried only hoop stress. Through an oversimplification of the 
geometry the erroneous conclusion was drawn that the tension in 
the cords would be constant from skirt to vent. It is interesting 
to note that in the analyses of Taylor the cloth stresses were in 
ther meridional direction only, while Stevens and Johns obtained 
the same shape by assuming cloth stresses to be in the hoop 
direction only. 
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Beck ^ summarized the German work on parachutes in 1942. Ex- 
pressions are given for membrane stresses in surfaces of revolution 
for constant and linearly varying pressures. 

Investigations of conical ribbon parachutes were conducted by 
Jaeger, Culver and Della-Vedowa of Lockheed Aircraft5-8 in 
1952. A method was developed for estimating the maximum hori- 
zontal ana radial tape loads in conical ribbon parachutes loaded 
by uniform pressure. The upper portion of the canopy was assumed 
to be a cone with zero meridional load and the lower portion was 
assumed to be an ellipsoid, with zero hoop tension. 

A comprehensive report by Topping, Marketos, and Costakos^ in 
1955 brought together all previous work in the field of parachute 
stress analysis. This work approaches the practical limit of 
refinement in the approach used, i.e., the fitting of functional 
curves to the parachute shapes. Analytical solutions are given 
for fully inflated shapes of solid flat, extended skirt, personnel 
guide surface and conical rlngslot parachutes loaded by uniform 
pressure. Various modifications to the Taylor equation were made 
to better approximate shapes observed in photograph of parachutes 
in steady descent. Radial cord tensions and cloth stresses con- 
sistent with the analytically determined shapes were calculated. 
Several representative dimensions obtained from the analytical 
solution for each of the four types of parachutes studied agreed 
with corresponding dimensions obtained from photographs within 
9 percent. For the flat canopy, dimensions agreed within about 
4 percent. Since this degree of accuracy might appear to be ac- 
ceptable, it should be pointed out that the variation in stress 
implied is much greater than the percentages given. For a typical 
nylon material (MIL-C-7020 Type I, cloth) 5 percent variation in 
strain corresponds to a 50 percent change in stress. The effects 
of varying such parameters as the suspension line length, number 
of gores and cone angles, as well as the effect of other details, 
were studied. 


6 


NVR-6432 



NOHTHItOP 

A generalized equation for Taylor type curves was given by 
Lester^ in 1962, 

sin 0 = \ - 1 + J 

By varying the parameter X, a family of curves is generated 
which may be used to represent different parachute types. 

11-14 

Heinrich, together with Monson and Jamison developed a 

method for computing canopy cloth stresses for a known profile 
and pressure distribution. Since the gore centerline is usually 
the profile observed in photographs, equations were developed 
using this profile. A set of five simultaneous equations are 
given which can tfe solved by an iterative process to obtain cloth 
stresses. This method is also presented in the "Parachute 
Handbook"^. 

1 f 

A dynamic analysis presented by Asfour in 1966 is based on the 
assumption that the major (circumferential) stress is a result of 
a dynamic force in a plane perpendicular to that of the opening 
shock. Circumferential stresses are computed by equating the 
energy absorbed by the cloth to the kinetic energy of the air mass 
expanding to fill the canopy during inflation. It is assumed 
that each horizontal member forms a circular ring around the 
canopy, i.e., the bulges disappear, at the time of maximum stress. 
This assumption, together with that of a linear stress-strain 
relationship, gives these rings linear spring rates and makes it 
possible to solve the energy equation in closed form. 

In the Apollo Block I stress analysis^ and '*' 0 Speakman and Topp 
treated the parachutes as surfaces of revolution. The pilot 
chute and the reefed drogue and main parachutes were assumed to be 
hemispheres loaded by uniform pressure. Analysis of the drogue 
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and main parachutes in the full open condition were improved by 
approximating the shape observed in drop test photographs by 
ellipsoids. In the analysis of the main canopy, further refine- 
ment was made by using a nonuniform pressure distribution based 
on steady state wind tunnel test results. 

In the Apollo Block II stress analysis of the main parachute by 

1 8 

Ranes, Topp and Utzman, a new method was used which accounted 
for th( effect of the bulging of the sails between the radial 
tapes and for the relationship between meridional curvature and 
forces. Pressure distribution data from wind tunnel tests were 
incorporated and the curvature of the radial tapes was assumed 
to be elliptical. Sail strain compatibility was used to select 
a best-fit ellipse. While this method represented a major step 
in development of a rational analytical method, it was limited 
by its reliance on a functional curve to describe canopy curvature. 
In effect, the canopy was forced into a convenient shape without 
regard for local member stiffness. As pointed out in that analysis, 
a change of 4 percent in the diameter at a particular section 
could result in a 50 percent change in the calculated stress. 

In the Block II (H) Program, a major step was taken in the de- 
velopment of a more realistic method for determining the load 
distribution in a parachute canopy. This method, which is given 
in Section 3, determines the unique shape and internal load distri- 
bution which satisfies equilibrium and boundary conditions for a 
given parachute and applied loading. The final stress analysis 
for the Apollo Block II (H)^9 used this method. However, verifi- 
cation of the method by correlation with test results was limited. 
Verification and continued refinement of the method by correlation 
with the test data obtained in the Apollo ELS test program is one 
of the goals of this study. 
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SECTION 2.0 

DISCUSSION OF PARACHUTE STRUCTURAL ANALYSIS 


The major factors encompassed in the structural analysis of 
a parachute are those appearing in the margin of safety equation. 


Ms = Ultimate strength x design factors ^ 
Limit load x safety factor 


Each of these factors is discussed in the paragraphs below. 


2.1 MARGIN OF SAFETY 

The margin of safety of a component is the ratio of excess 
strength to the required strength of the component. An ideal 
design would have a zero or small positive margin of safety 
for every member. 


2.2 SAFETY FACTOR 

The safety factor is the ratio of the required ultimate strength 
to the limit load applied to a member under the worst case that 
can be expected to occur with reasonable probability within the 
envelope of normal operating conditions. This is the "ignorance" 
factor which covers the uncertainties in load and strength pre- 
dictions. Since the weight and volume of the parachute will vary 
almost linearly with this factor, it is desirable to use the 
lowest factor consistent with the accuracy of the analyses. 

V 

2.3 DESIGN FACTORS 

Design factors account for reduction in ultimate strength due 
to environmental degradation and for increases in some member 
loads due to nonsymmetrical loading of the parachute. 
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2.4 ULTIMATE STRENGTH 


The ultimate strength cf a parachute member is usually limited 
by the strength cf splices within the member or at Joints with 
other members. In the Apollo program, representative joints 
were tested to destruction in the laboratory to determine the 
expected strength of the joints relative to the base material 
strength. 

Joint efficiency = Joint strength 

Base material strength 


and 


Ultimate Strength = Material Rated Strength x Joint Efficiency 


Strengths based on static (low speed) joint and material tests 
are presently used for most members. For pilot parachute risers, 
however, there is evidence that the high load onset rate at 
time of canopy stretch gives important dynamic effects. Section 
5 is a study of the behavior of the Apollo pilot chute riser 
under dynamic loading. 


2.5 LIMIT LOAD 

The limit load in a member is a function of the external loading 
on the parachute and of the internal distribution of this loading. 


2.5.1 External Loading 

Primary external loads on a parachute are the riser load and 
the reacting differential pressure on the canopy surface. Most 
of the effort in parachute loads analysis is directed toward pre- 
diction of riser loads, and methods have been developed which 
predicts this load with a high degree of accuracy. (See, for 
example. Volume 1). 

There is very little quantitative information on the differential 
pressure distribution on parachute canopies. The most meaningful 
empirical data available comes from the work done by Melzig. 
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References 20 and 21 give the results of wind tunnel and free 
flight tests in which local pressures were measured at four 
points in parachute canopies during opening. Further dis- 
cussion of the details of these tests is given in Section 6. 

Figure 1 is a plot of the data for a typical wind tunnel run 
from Reference 20. The variation in differential pressure coef- 
ficient versus time at the four pressure transducers is shown. 

It is seen that a pressure peak starts at the vent and moves 
.’apidly down the canopy to the skirt. Correlation with load- 
time traces show that maximum load occurs at approximately 
the time this peak reaches the skirt, although development of 
maximum canopy projected area occurs somewhat later. 

The results of drop tests given in Reference 21 differ from the 
wind tunnel data in that the pressure peak never reaches the 
skirt because of the decay in the dynamic pressure during opening 
under finite mass conditions. Canopy loading was very low (1/3 psf) 
for the tests of Reference 21. It would be expected that behavior 
of the more highly loaded spacecraft parachutes would be some- 
where between the limits explored in References 20 and 21. 

For reefed parachutes there are no pressure measurement data 
available. In the stress-time study of Section 4.3.4, the feasi- 
bility of Inferring the pressure distribution from the shape of 
a parachute is explored. 

2.5.2 Internal Loading 

Loads in the Individual parachute members for a given riser 
load and pressure distribution are determined by the method given 
in the following section. 
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Differential Pressure Coefficient 
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SECTION 3 

PARAC -JTE INTERNAL ID AD DISTRIBUTION 

The determination of the Internal load distribution for a given 
parachute and set of external loads is the subject of this section. 

The basic analysis method, which has been implemented in a com- 
puter program called CANO, Is derived in Section 3-1 Ihis 
analysis is used for most of the test data correlation work In 
Section 4. A further refinement of the analysis, which is im- 
plemented In Program CANO 1, Is presented in Section 3-2. 

3.1 PROGRAM CANO - AN INTERNAL LOAD ANALYSIS FOR PARACHUTES 

The analysis method which follows determines the unique shape 
and Internal load distribution which satisfies equilibrium and 
boundary conditions for a given parachute under the influence 
of known riser and aerodynamic forces. 

This analysis it applicable to either reefed or fully open 
symmetrical designs, such as the ribbon, ringslot, and ringsail 
parachutes, which have meridional members which may be assumed 
to carry all meridional forces. 

A general description of the method and assumptions is given 
first, followed by derivation of the equations and a solution 
algorithm. 

3.1.1 Description of the Method 

The parachute is treated as a deformable membrane, and finite 
elements are used to represent the structure. The structural 
model of the canopy consists of horizontal slices (or segments) 
conveniently spaced alor^ the meridional members from the skirt 
to the vent. Each segment is defined by its length along the 
meridian, the free length o" the horizontal member connecting 
adjacent meridional members, and the load strain curves for these 
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members. Loads and geometry are computer at reference stations 
located at the midpoints of the segments. 

The following assumptions are made: 

1) Meridional curvature is constant over each segment, 
and the curvatures of adjacent segments are tangent 
at their junctions. 

2) All meridional forces are In the meridional members 
vradjal tapes). 

3) Horizontal members have no meridional curvature, but 
each has a constant radius of curvature which Is 
normal to the meridional member at the reference 
station. 

The simplifying assumptions given above make It possible to 
compute the end reaction of a horizontal member as a free body. 
This end reaction is then resolved into three mutually perpendi 
cular forces tangent and normal to the meridional member, from 
which the loading and radius of curvature of the meridional 
member can be computed. 

The variation of differential pressure along the meridian is 
represented by a nondlmensional distribution curve. Magnitude 
of the differential pressure is adjusted to give overall force 
balance by two factors; 1) the average pressure acting on the 
projected area of the skirt, and 2) an Integration factor 
which accounts for the shape of the distribution curve and the 
geometric porosity of the model . 

A solution is obtained for a trial skirt diameter and pressure 
integration factor by establishing equilibrium geometry, first 
at the skirt then at each successive station to the vent. This 
trial solution is then examined for compatibility with boundary 
conditions at the vent. (See discussion below) This process 
is repeated with corrected trial skirt diameters and pressure 
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factors until the boundary conditions at the vent are satisfied 
and an overall force balance is achieved. 

Boundary conditions at the vent are checked as follows: The 

load in a meridional member at the vent band station is reacted 
by a vent line. The stretched ? ^th of the vent line under 
this loading is compared with tl - computed diameter. The other 
boundary condition is zero slope of the meridian at the vent, 
if an idealized case in which no pressure is acting on the vent 
lines is considered. This singularity is avoided in a practical 
manner by allowing a small positive slooe. This slope is limited 
by a requirement that the total vertical component of the meri- 
dional member loads not exceed 75 percent of the force .hat would 
be produced by the local pressure acting on an area equal to 
that of the vent. 

3.1.2 Derivation of Equations 

Since the solution procedure is not straightforward, no attempt 
is made to follow the computing sequence in the derivations 
which follow. A solution algorithm is shown in the form of a 
flow diagram of computer program CANO in Figure 2, page 26. 
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Equilibrium below the skirt - unreefed 

A profile view of a fully inflated, unreefed parachute Is shown 
below. 



For a given skirt diameter, equilibrium geometry is established 
as follows: 

a = sin’^/Lgfl + ^)] (1) 


L M cos a 

where: a = suspension line convergence angle 

L' = initial length of suspension lines 
s 

* 1 , = strain of the suspension lines under load, Pl 
V = applied riser load 

M = number of suspension lines (also number of gores) 

Equations (l) and (2) are solved simultaneously by the method 
shown in Figure 3, Page 27. 
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Equilibrium below the effective skirt, reefed 

For reefed parachutes It Is assumed that differential pressure 
acts on the Inflated portion of the canopy only, so that the 
radial tapes are straight In the unlnflated portion. For a 
given development angle, b, equilibrium geometry Is established 
as follows: 


I 



Equations (3) through (8) are resolved simultaneously by the 
method shown in Figure 4, Page 28. 
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From View A -A 




(9) 

( 10 ) 

( 11 ) 


A = M(Ag + A^) sin 0 


( 12 ) 


Local differential pressure 



(13) 


where: p = local differential pressure 

p' = pressure coefficient from pressure r ’ ^tribution 
curve 


kp « Integration factor-a function of the shape of the 
distribution curve and geometric porosity 

The above equations are solved simultaneously by the method 
shown in Figure 3 or 4. 
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EQUILIBRIUM OF A TYPICAL HORIZONTAL ELEMENT 


R 



View In Plane A 

Prom summation of forces in direction R' 




N* = -E& 

2 sin ^ 

(14) 

where : 

N^, = unit 

tension load in horizontal member 




h = 2r sin tt/m 

(15) 

Therefore 

N* - P- r - s l»- TT / M 

sin t 

(16) 

The arc 

length 3 l 

is given by 




S L * = 2 ♦ Rs 

(17) 


Since: 


Rg sin t = r sin tt/M (l8) 

S L '= 2^ r sin tt/M (19) 

sin ♦ 


The length of the horizontal element 

S L = S M (1 + e s ) (20) 

where: Sm = manufactured length of the horizontal 

member 

e s = unit strain of the horizontal member for 
° load Nty 

A value for ^ is found by iteration to give 

S L = (21) 

This procedure is shown schematically in Figure 2, page 26. 
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Resolution of Niji 


The edge force In a horizontal member, , Is resolved into 
three mutually perpendicular components; (l) APp, tangent to 
the meridional member, (2) Np, normal to the meridional member 
in a plane which includes the central axis, and (3) , a cir- 

cumferential force. 

From the view in Plane A: 

Nq = N,^, sin 'i'/cos u (22) 

and N p = N.,. [cos 'i + sin sin a | (23) 

r y cos a 

Angle a is defined as follows 
From Equation (15) 
h = 2r sin n/M 

sin ci = - = — s - n - = (sin n/M) sin 0 (24) 

2Rq 2r/sin 0 

Np is resolved into components in the desired directions by 
first resolving it into two components in the horizontal plane 
as shown below. 



View In Horizontal Plane 

Nq = Np cos n/M = Ny cos n/M [cos + 3 ^ n ~ s A n - (25) 

cos a 

Nrh - Np sin n/M (26) 
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The resolution of Nr r into the directions of N R and Ap R 
be seen with a view in the plane of a meridional member. 


can 



From u F z = 0 


View In Plane Of Radial Tape 


APp = 


jjg_» cos 0 


sin 0 


(27) 


From ^ Fr = 0 and equation 23 and 24 

N r * = Nr h sin 0 = Np sin a 

= Nt sin a L cos i + sin t sin a/cos a J 


(28) 


Combining equations (22) and (28) 

N R = N Q - Nr* 

Nr= N^sin */cos a - sin a (cos * + sin <' sin a/cos a)J (29) 
Substituting equation (28) into (27) 


APr = Nij- sin tt/m cos 0 (cos ^ + sin ♦ sin a/cos a) 

(30) 

■0 
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EQUILIBRIUM OF A SEGMENT OF A MERIDIONAL MEMBER 
(loaded by two adjacent horizontal member ends) 



A0 is sufficiently small to make APr AS negligible compared to 
P R . Summation of forces in direction R results in: 

2(2 N r R 0 sin A0/2) = 2 (P r sin A0/2) 

therefore 

R0 = P r /2N r (31) 


The above equations can be combined to give a form similar to 

/ N0 NO \ 

the general membrane equation (p = r^ + ^ ): 


P 


-[ 


R 


N 


0 sin 


2r R* sin 1 
v M 


r cos 




M 


sin 2 — sin 2 0 
M 


/ 
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LOAD IN THE MERIDIONAL MEMBER, Pr 


At the skirt the axial load in the meridional member is equal 
to the suspension line load. 


P R (skirt) - - 

K M sin 0 


(32) 


At subsequent stations P R is computed by subtrac' lng the ac- 
cumulated AP^ from the initial value. Since there are two 
horizontal member ends acting on each meridional member, the 
load at station j is given by 


p Rj ' P R (J . 1) - < 2 ip R AP R AS/2 )j <33) 

Length of segment J is given 

ASj = AS» (l+e Rj ) 

where: t 

AS j is the manufactured length of the segment 

e„ is unit strain under load Pr,. 

J J 
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3.1.3 Solution Algorithm 

A flow diagram for Program CANO is presented In Figure 2. 

A user's manual which includes a listing of this program is 
given in Reference 22. 


25 


NVR-6432 



v_c-<~"T3; ry , strain 
Pressure : irve. R 


ASSUME SKIRT l'IA, r , cr LEVEL'. : 


ASSUME PRESSURE COEFi 


Reef ea 


RT E Q UILIBR IUM! IGKIRT ESUIJ 


assume mer]?:u;ia:, radius 

OF CUPVATUR. . R 


rtSSUIiE SAIL BULGE ANGLE, 4 


COMPUTE SEGMENT GEOMETRY AND LuADI 


ICOMPARE HORIZONTAL RlBBt N LENGTHS 


CHECK V FOi DIVER GEN C 


COMPARE h > 


COMPARE Pu 


CHECK STATION NO., : 


i=N 


CHECK AXIAL LOAD BALANCE 


K VENT DIAMETER 


IPRjNT RESULTS I 


I low Diagram for Program CANO 


NVR 04 3 P 





















Fig. 3 . Flow Diagram Detail Showing Skirt 

Equilibrium for an Unreefed Parachute 
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ASSUME TRIAL BULGE AREA, A 


COMPUTE LOCAL PRESSURE 


ASSUME TRIAL BULGE ANGLE 


COMPUTE HORIZONTAL LOAD, N*, AND LENGTH, S 


, I 

COMPUTE ARC LENGTH, Sj 

, I 

COMPARE S L AND S £ 

, r 

COMPUTE BULGE AREA, An 

- -I , 

COMPARE A R AND An 

, 3 

COMPUTE J0, A0, P R , AP R , e, d, g, h 


Fig. 4 . Flow Diagram Detail Showing Skirt 
Equilibrium for a Reefed Parachute 
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3.2 PROGRAM CANO 1 - AN INTERNAL LOADS ANALYSIS FOR RIBBON 
PARACHUTES 

In deriving the equations for canopy Internal load distribution. 
Section 3-1* it was assumed that the horizontal members bulge 
outward perpendicular to the meridional members, i.e., angle 
0=0 in Figure 5 . It is seen in the photograph in Figure 

5 that this assumption is not realistic for ribbon parachutes. 
The vertical members between the meridional members, or radial 
tapes, pull the skirt band and adjacent horizontal members up so 
they are not perpendicular to the meridians. As a result, the 
vertical members pick up axial load which would otherwise be in 
the radial tapes. This analysis attempts to account for the 
effect of the vertical members on the internal load distribution 
in the canopy. 

3.2.1 Description of the Analysis Method and Assumptions 

This analysis is identical to the CANO analysis given in Section 
3.1, except that the equations and solution procedure are modi- 
fied to include the effect of the vertical members. The struc- 
tural model consists of horizontal slices (or segments) of the 
canopy. Each segment is represented by a length of radial tape, 
a horizontal member, and a vertical member which is uniformly 
distributed across the width of the gore. Loads and geometry 
are computed at reference stations located at the midpoints of 
the segments. 

The following assumptions are made: 

1) Meridional curvature is constant over each 
segment, and the curvatures are tangent at 
the junction between adjacent segments. 

2) The horizontal ribbons have no meridional 
curvature. 

3) The edges of a horizontal ribbon lie in planes 
parallel to the Plane B (Figure 6. ) 
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HORIZONTAL RIBBON 


VERTICAL RIBBON 


SUSPENSION LINES 


Figure 5 • Ribbon Parachute With Vertical Members 
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4) The projection of the horizontal ribbon on 
Plane A (Figure 6) is a circular arc. 

5) The horizontal ribbon is held in equilibrium 
by a variable distributed force from the; 
vertical ribbons. This force is parallel to 
the horizontal ribbon width, and is of a 
magnitude such that the resultant of the force 
system acting on the horizontal ribbon is in 
Plane B (Figure 6). 

6) The axial loads in the vertical ribbons are 
transferred to the radial tapes over each seg- 
ment from skirt to vent in proportion to the 
decrease in gore width. 

3.2.2 Derivation of Equations 

Since the solution procedure is not straightforward, no 
attempt is made to follow the computing sequence in the 
derivations which follow. A solution algorithm is shown in 
the form of a flow diagram for computer program CANO 1 in 
Figure 7 , page 41 . 


31 


NVR-6432 







EQUILIBRIUM OF A TYPICAL HORIZONTAL RIBBON 


From the assumptions that: 

1) The load increment ANy (from the vertical 
ribbon) is perpendicular to Plane A. 

2) The pressure load is in Plane A. 

3) The resultant of the two loads is in Plane B. 

It follows that 


ANy = p cos § tan B (36) 



N S 

N T 


N 


T 

= ^ p cos § Rg d§ = p Rg sin $ 
o 

«= \ p cos § tan B Rg d § = p Rg sin 'll tan B 


¥ 

P = pR S t 1 " ^ sin 5 d§ J ® P R g cos 

A° 


N = [Ng 2 +N T 2 + (Np' “ pr sin tt/M [esc 2 ^ + tan 2 b] 


■\i 


(37) 

(38) 
(39) 
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Length of the elliptical arc Is given by 

S' = V(1 - a in 2 3 sin 2 ?) d § 

^ cos p i 

* 

n/ 2 

, = 2 r sin tt/m ^ 


St' ■ 


sin t cos P V 

V 


sin 2 9 sin 2 = ) 


( 40 ) 


Stretched length of the horizontal ribbon is given by: 


Sl - S M ( 1 + e s > 


(U 1) 


where : 

Sj^ = manufactured length of horizontal ribbon 
= unit strain for load 

A value for $ is found by iteration to give 

S L - S £ (Its) 

This procedure is shown schematically in Figure 7, page 41. 

The three components of the horizontal ribbon end reaction, 

Ngi Np 1 , will be resolved individually into three orthogc. 
components perpendicular and tangent to the radial tape. 
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RESOLUTION OF REACTION Nt 


i 



View in gore 
centerline plane 



View in horizontal plane 


i 



forces on radial 
from one horizontal 

View in radial plane 


sin 0 ' - r 008 n/M 

J (r cos Tt/M) 2 +(r/tan 0) 2 

_ sin 0 cos n/M 

(43) 

cos a 


cos 0’ * Jl - sin 2 0' = 

cos 

_0 

a 

(44) 

a = sin”* (sin 0 sin n/M) 



Ntjj = N t cos 0' 

(45) 

NTy = N(j sin 0' 

(46) 

N©g - N^ijj sin n/M 

(47) 


NRr “ Nt V cos 0-th 008 n / M sin 0 

Nr B - 0 (48) 

AFRg * N»]ty Sin 0+Nfu 008 n /M C 08 0 

(49) 
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Resolution of Np and Ng 



i 

View in plane A 


Resolution of Reaction Np 



View in horizontal plane 

N© a = Np cos tt/m (52) 

Nr h = Np sin n/M 
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(53) 



I 

(Forces on radial tape 
from one horizontal ribbon) 


View in plane of radial tape 


Combining components of Nt, Ns, and Np’ 

N 9 " % + N ©B 

From Equations (37), (38), (39), (44), (45), (47), (50) and (52) 

N q = p . r s ^ n .- i rr / . M [ cos tt/m (cos t + sin $ tan a) 

9 sin ♦ L 


. tan » cos 0 J 

cos a 


(55) 


n r = n Ra + Nr b 

From Equations (37), (39), (48), (50), (52), and (53) 


Np . pr sln - n/M j~ Bln -♦ 


- sin a (cos ♦ + sin h tan a)j (5 6) 


sin 'll 1 cos a 

AP r * 4p Ra + APRb 

From Equations (37), (38), (39), (43), (44), (45), (46), (49), 
(50), (52), and (54) 

APr « pr si fl , n / M i S i n tt/m cos 0 (cos ♦ + sin ♦ tan a ) 
sin V *- 

+ tan 9 cos tt/m ■ 1 (57) 

ccs a J 
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The load In the vertical ribbons is zero at the skirt since this 
is a free edge. Each horizontal ribbon adds an increment of load, 

MyAg, so that the loading at station J is given by 

Nv J = + (A% AS/2) (J-1 ) + (AN V AS/2)j (58) 

Since the verticals intersect the radials as the gore width de- 

creases from skirt to vent, the vertical ribbon load is transferred 
to the radial tape. It is assumed that this transfer is proportional 
to the decrease in gore width between reference stations. The force 
applied to a radial over a length AS' is given by: 



&P Rv AS' = N v AS M cos tt/M 

APr„ = N„ — — cos tt/m 
v v 4 S' 


^ (approx) 

Flat Pattern of Gore 


(59) 
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The load In a radial tape at station (J ) Is given by: 



p R (j) - p R(>i) - 2 < 4P r As/2 >(j-i) - 2 < ap r as/2 >j 


( 60 ) 


Curvature of Radial Tape 

The following approximation Is used for the radius of curvature 
of the radial tape: 

R .ajjvfs (61) 

© 2 % 

This equation is identical to that derived in Section 3.1 except 
that the term NvSm is added to account for the load in the verti- 
cal members. In the solutions obtained for the drogue chute 
(Section 4.1.6), use of this equation results in equilibrium for 
a vertical fcection through the canopy when checked by the method 
shown below. 
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0 = 



( 62 ) 


Solution Algorithm 

The equations for skirt equilibrium, radial tape geometry and 
differential pressure are the same as those given in Section 3.1. 
A solution algorithm is given in the form of a flow diagram for 
computer program CANO 1 is given in Figure 7. 


40 


NVR-6432 




41 


NVR 6432 













SECTION 4.0 


TEST DATA REVIEW AND APPLICATION OF INTERNAL LOADS ANALYSIS 

The objective of this section is to upgrade the structural anal- 
ysis methods by use of the test data acquired during the Apollo 
Test Programs. The forms of test data available are: l) high- 

speed motion picture coverage from onboard, air-to-air, and 
ground cameras; 2) load -time traces at the parachute riser; 

3) dynamic pressure history from ASKANIA tracking; and 4) failure 
analyses of drop tests In which parachutes were damaged or de- 
stroyed. Appendix A summarizes the drop tests of the Apollo ELS 
Block I, Block II and Block II (H) programs. 

In this section, two methods are used to verify the analysis 
methods presented In Section 3* First, drop test photographs 
are compared with the shape predicted analytically for the test 
loading conditions. Since a unique shape is predicted for a 
given pressure distribution and axial loading for a given para- 
chute construction, correlation of the predicted and actual 
shapes is evidence that the internal load prediction is correct 
(if the pressure distribution used is correct). This is further 
verified by comparison of predicted ultimate strength and failure 
modes with test results. 

Additional analyses are presented to show the sensitivity of 
canopy size and internal loading to variations In canopy con- 
struction and pressure distribution. 

Analysis of the main parachute is extended to include a stress- 
time study throughout the opening process and optimum weight 
calculation for the spacecraft design. 

Three types of parachutes are used in the Apollo Earth Landing 
System, the drogue, pilot, and main parachutes. The first 
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two parachutes to be deployed are drogue chutes, which are ribbon 
type parachutes that incorporate one stage of reefing. Three 
pilot parachutes are deployed just as the drogues are discon- 
nected to extract the main parachutes from their packing bags. 

The pilot chute has a ringslot type canopy with an overinflation 
line attached to the skirt. Final deceleration of the space- 
craft to descent velocity is accomplished by the three main 
parachutes. These ?.re ringsail parachutes with two stages of 
reefing. The analyses that follow will consider the parachutes 
in the order of deployment. 
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4.1 DROGUE CHUTE 

4.1.1 Drogue Chute Structural Model 

Structural details of the drogue chute are given in Figure 8 
and the corresponding structural model is defined in Figure 9. 

In modeling the canopy, all horizontal members, including the 
skirt and vent bands, are represented as finite elements. Initial 
ler^ths and tape rated strengths are taken from the engineering 
drawings and appropriate load-strain curves are input to the 
computer in tabular form. The reefing and overinflation lines, 
which are attached at or near the skirt band, are Included in 
the model by the method discussed in Section 3» 

Modeling of the meridional members presents a problem because 
of vertical ribbons between the radial tapes. There is experi- 
mental evidence that the verticals carry significant meridional 
load, but the method used in program CANO assumes that all the 
meridional load Is carried by the radial tapes. To approximate 
their effect in the structural model used by CANO, the stiffness 
contributed by the verticals is added to that of the radial 
tapes. Results of a study to determine the sensitivity of the 
solution of a fully inflated canopy to variations in meridional 
stiffness are given in Figures 10 , 11 £.nd 12 . Canopy profile, 
horizontal member loading# and radial tape loading for three 
models differing only in the stiffness of the meridional members 
c.re compared. Model I has radial tapes only; Model II has the 
two verticals nearest the radial fully effective; and Model III 
has all of the verticals fully effective. In the latter two 
models, the vertical stiffness is diminished linearly from fully 
effective at the station with maximum diameter to zero at the 
skirt. It is concluded from these figures that modeling variations 
of the meridional members are not significant in determination 
of the canopy profile but can affect the maximum load in hori- 
zontal tapes. A difference of approximately eight percent in 
the maximum horizontal tape load Is observed between Models I 
and III. Further study of the effect of vertical ribbons using 
the analysis of Section 3-2 is presented in Section 4-. 1.6. 
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Pig. 9. . Drogue Chute Structural Model 
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Fig. 10 iiroguf Chute Canopy Profiles Ter Various Structural Models 
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4.1.2 Drogue Chute Pressure Distributions 

The drogue parachute Is analyzed for both the reefed and open 
stages since not all canopy members are subjected to their 
critical loads in the same stage. 

Since pressure distribution data for reefed parachutes is not 
available, the canopy is assumed to have a uniform pressure 
differential over the inflated portion of the canopy'. The 
degree of canopy inflation depends on the length of the reefing 
line, which restricts the inflow of air, and on the porosity of 
the canopy. The degree of inflation of the reefed drogue was 
determined from photographs taken during an aerial drop test. 

(See Section 4.1.4) Figures 13 and 14 show the horizontal 
ribbon loading, and radial tape loading for the reefed drogue 
chute with a uniform pressure distribution. To determine whether 
this is the appropriate pressure distribution, the predicted 
canopy shape is compared to the measured canopy shape in Section 
4.1.4. The reefed canopy is analyzed for a riser load of 21,000 
pounds which is the peak canopy loading occurring in the flight 
test used for the shape comparison. 

Some empirical data is available for pressure distributions in 
ribbon canopies for the full open condition. The wind tunnel 
pressure distribution data. Reference 20, is applicable for the 
drogue chute since the dynamic pressure decay is only about 
2 percent from time of dlsreef to maximum load. Test data for 
a circular ribbon parachute with 18 percent geometrical porosity 
are given in Reference 20. The drogue chute is a conical ribbon 
parachute with a geometrical porosity of 22 percent. The drogue 
chute has an overinflation line, which might affect the pressure 
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rt Station Number 

Fig. 13. Horizontal Ribbon Loading for the Reefed Drogue Chute 
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Fig. 14. Radial Tape Loads for the Reefed Drogue Chute 


muinmf 

distribution at the skirt, and it Inflates from a 43 percent 
reefed condition, whereas the wind tunnel model Inflated from a 
zero percent reefed condition. Mo test data are available to in- 
dicate the effects of these differences. 

In Figure 15* the data of Reference 20 is replotted to show the 
differential pressure across the canopy at different times for 
the four test velocities. Two significant observations can be 
made from Figure 15; 1) the general shapes of the curves are 

similar at the three event times shown--maximum diameter, maxi- 
mum load, and steady state, and 2) no trend with varying velo- 
city is indicated. On the basis of these observations it is 
concluded that the mean curve indicated in Figure 15 is a reason- 
able approximation for use in the analyses which follow. 


In order to determine the effect of variations in pressure distri- 
bution on canopy shape and internal loading, a sensitivity study 
was performed using the four pressure distributJ ns shown in 
Figure 16 . Distribution A is the mean curve from Figure 15 and 
curves B, C, and D are arbitrary distributions. In distribution 

B, the pressure is constant throuhout the canopy. In distribution 

C, the maximum pressure is at the vent while in distribution D 
the maximum is at the skirt. In distributions C and D the maxi- 
mum pressure is twice the minimum value. 

Figure 17 shows the variations In the canopy profile for the 
four pressure distributions. As expected, the peak pressure 
at the skirt gives the largest canopy diameter, while peak pres- 
sure at the vent gives the largest profile height. 

The variation of horizontal ribbon load with these pressure 
distributions is given in Figure 18 . It Is seen that distri- 
bution D (peak pressure at the skirt) gives peak loading 11 percent 
higher than that of distribution C (peak pressure at skirt). The 
peak horizontal tape load occurs higher in the canopy. as the peak 
pressure moves from skirt to vent. 
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Fig. l'j , Pressure Distributions At Various 
Time 3 After Deployment 
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Fig. 16. Canopy Pressure Distributions 
Used for Sensitivity Study 
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Station Number vent. 






In Figure 19 the radial tape loads for the four pressure 
distributions are shown. Distribution D produces the highest 
loading at the skirt, while distribution C gives higher loads 
in the crown. 

It is concluded from these figures that the differential pressure 
distribution across a parachute canopy has a first order effect 
on internal loading. Further investigation of pressure distri- 
butions in inflating canopies is, therefore, essential for ac- 
curate internal load predictions. 
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Fig. 19 . Radial Tape Loads in Drogue Chute 
for Various Fr»- isurt- Distrioutions 












4.1.3 Drogue Canopy Growth 

Results of a study of the variation in canopy shape and internal 
loading with applied axial load are shown in Figure 20 . So- 
lutions were obtained for the disreefcd drogue chute, structural 
Model II, for axial loads of 5*000, 10,000, 15,000 and 23,340 
pounds, with all other parameters constant. In Figure 20 
maximum diameter, projected area, radial tape load, and hori- 
zontal ribbon load are plotted versus applied axial load. The 
variati n of these quantities is practically linear with load. 
The slight curvature observed can be attributed to the shapes 
of the load strain curves of the parachute materials. Since 
this effect is negligible for small changes in load, canopy 
loads can be proportioned directly with axial load. 
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Fig. 20. Variation of Canopy Size and 

Internal Loads with Axial Load 



4.1.4 Drogue Chute Film Analysis 

Three types of photographic coverages are available from the 
Apollo Drop Test Programs: air-to-air, ground-to-air, and 

onboard. It is desirable to measure the profile view of the 
canopy in order to obtain the most accurate shape comparison. 

These measurements would have to come from the air-to-air or 
ground-to lir coverages. Due to the size, velocity, and de- 
ployment altitude of the drogue chute, the air-to-air and ground- 
to-air coverages are inadequate to provide shape measurements. 

It is, therefore necessary to use onboard cameras coverage. 

The individual ribbons can be identified in these photographs 
so that canopy diameter as a function of station number (as 
defined in Section 4.1.1) can be measured. To obtain actual 
dimensions of the canopy from the photographic images, the 
measurements are converted as follows: 

Diameter = Image Diameter x Lens Factor x Distance from Camera 

where the lens factor is calculated by measuring the image of 
a known length at a known distance. The distance from the 
camera to the skirt is calculated from the free lengths of the 
risers and suspension lines and their elongations. Distances 
from the skirt to specific horizontal members are taken from 
the analytically predicted geometry. 

Figure 21 shows the comparison of the reefed drogue shapes 
for constant pressure over the inflated portion of the canopy. 

From the onboard photographs, it was seen that the first eleven 
horizontal tapes above the skirt were not inflated. The pressure 
distribution used to predict the analytical 3hape was, therefore, 
adjusted so that the first eleven horizontals had zero differential 
pressure. Since the predicted canopy shape is close to the ob- 
served shape, the computed internal canopy loading may be assumed 
to be close to actual values. (See Figures 13 and 14 . ) 
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Station Number 

Comparison of Predicted and Measured 
Reefed Canopy Shapes for the Drogue Chute 
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Figure 22 shows the comparison of the disreefed drogue shapes 
for pressure distribution A in Figure 1 6 . The effect of the 
overinflation line on the disreefed canopy is to control the 
expansion of the skirt and not to restrict the inflation of the 
canopy. All of the horizontal tapes, therefore, become inflated. 
It is seen in Figure 22 that good agreement exists in the upper 
portion of the parachut. but in the lower portion the predicted 
diameter is too large. A significant difference in shape exists 
at the skirt which would result in larger than actual loads in 
the radials and overinflation line. This discrepancy is attri- 
buted to the presence of the overinflation line which causes a 
large development angle at the skirt (25°). This would tend to 
reduce the differential pressure in the skirt area, so that the 
test data used to obtain this distribution curve is not appli- 
cable. 
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4.1.5 Failure Analysis 

Table 1 gives a comparison of predicted strengths and failure 
modes with test results for the two tests In which ultimate load 
was applied to fully opened drogue canopies. This comparison, 
although limited in scope by the number of test points available, 
shows reasonable agreement between predicted and actual strengths 
and failure modes. 

Table 1 Drogue Chute Failure Analysis 


Test 

No. 

Measured 

Riser 

Load 

Test Damage 
or 

Failure Mode 

Calculated 

Ultimate 

Strength 

Predicted 

Failure 

Mode 

84-4* 

27,250 lb 

No significant 
Damage 

30,300 lb* 

Gore 

Split 

** ^ 

99-2 

35,930 lb 

Gore split 
from Vent to 
Skirt 

, _ ** 
34,960 lb 

Gore 

Split 


* R8155-507 had 9 rows of stitching to secure the hori- 

zontals at the radials (Joint Efficiency = .75) 


R8155-505 had 10 rows of stitching to secure the hori- 
zontals at the radials (Joint Efficiency = . 83 ) 
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4.1.6 Analysis of the Drogue Chute Using CANO 1 

As was pointed out In Section 4.1.1, there Is experimental evidence 
that the vertical ribbons carry a significant portion of the meri- 
dional load in the canopy. To include the effect of the vertical 
members the CANO 1 Program was developed, as discussed in Section 
3.2. Solutions have been obtained for the reefed and fully open 
drogue parachutes to quantitatively evaluate the effect of the 
vertical members. The pressure distributions used are identical 
to those used to analyze the parachute by the CANO Program as 
discussed in Section 4.1.2 and 4.1.4, The riser loads have been 
revised to the design ultimate loads of 23,530 pounds for the 
reefed condition and 24,560 pounds for the fully inflated condition. 

Comparison of canopy profiles for the reefed and open parachutes 
are shown in Figures 23 and 24 , respectively. The main effect 
that the verticals have on the canopy shape is to lift the first 
few inflated horizontal ribbons so that they are not perpendicular 
to the radial tape as was assumed in the previous analysis. Radial 
tape profiles calculated for the reefed drogue canopy with and 
without verticals are nearly identical, as shown in Figure 23 . 

The gore centerline profiles differ only in the area of the first 
few inflated horizontal ribbons. Canopy profiles for the open 
drogue canopy also compare closely except in the area of the 
skirt. While the effect of the verticals on canopy shape does 
not appear too large, their effect on internal loading is sub- 
stantial. Figures 25 and 26 show the variation of hoop load 
in the horizontal ribbons and meridional load in the radial tapes 
for the drogue canopy with and without vertical ribbons. Figure 
25 shows that the hoop stresses rise sharply at the first in- 
flated ribbon in the reefed condition and at the skirt in the open 



NORTHROP 

condition due to the "hiking " of the skirt by the vertical 
ribbons. The contribution .he vertical ribbons in carrying 

meridional load is clearly . .:uicated in Figure 26 by the re- 
duction in radial tape loading for Y. th the reefed and open con- 
ditions. Figure 27 shows the v <cion over the canopy of 
the vertical ribbon membrane loading and the horizontal ribbon 
displacement angle for the reefed and open parachute. 

It should be pointed out that arbitrary assumptions were made 
in this analysis, so that supporting test data is necessary 
before it can be considered a design tool. 
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Canopy Radius, in. 

Fig . 24 . Disreefed Drogue Chute Profiles 
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Fig. 25 . Horizontal Ribbon Loading for Reefed and Open Drogue Canopies 
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skirt Station Number vent 

Fig. 26 . Radial Tape Loading fcr Reefed and Open Drogue Canopies 
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Fig. 27. Vertical Ribbon Loading and Horizontal Ribbon 

Displacement for Reefed and Open Drogue Canopies 











4.2 PILOT PARACHUTE 

4.2.1 Pilot Parachute Structural Model 


The pilot parachute has a twelve gore, ringslot canopy that is 
constructed as shewn in Figure 28 . A typical gore consists 
of five sails (42 lb/in cloth), three verticals (2-70 lb tapes), 
vent band and skirt band. The leading and trailing edges of the 
sails are reinforced by 0.6 inch wide seams. The radlals are 
reinforced by a 250-1 'd tape extending from the skirt to the 
trailing edge of Sail 4. The drag area of the canopy is con- 
trolled by an overinflation line which is attached to radial 
tape loops below the skirt. All meridional load in the canopy 
is assumed to be taken by the radial tapes and the effect of 
the vertical members is neglected. 

4.2.2 Pilot Parachute Pressure Distribution 

Pressure distribution data has been obtained from wind tunnel 
tests of ring3lot canopies by Melzig cv . The pressure distri- 
bution shown in Figure 29 was derived from this data. The 
predicted shape and internal loads for this design ultimate 
riser load of 4l85 lb and this pressure distribution were com- 
puted by means of program CANO. Figure 30 shows that the maxi- 
mum load to strength ratio for the sails occurs in the middle 
of Sail 3 at a value of 0.664. Also shown in tills figure is the 
variation in hoop load from skirt to vent. The hoop load Is 
largest at the leading and trailing edges of the sails. Due to 
the I creased strength in these areas (three plys of fabric in 
the leading and trailing edge seams) the ratio of hoop load to 
rated strength is lower than at the center of the sail Danel. 
Figure 31 shows the variation of radial tape load from skirt 
to vent. Maximum radial tape loading of 433 pounds occurs eight 
inches above the skirt. 


f 


74 


NVR-6432 



Vent Line 
425 lb Cord 

Vent Band 
2-900 lb Tapes 


Sta. 26 


Radial 
2 - 350 lb Tapes 


0.6 Seam 

(Typ.) 



250 lb 
Reinforcement 


2-70 lb 
Tapes 


-Overinflation Line 
650 lb Ccrd 


42 lb Cloth 
(Typ.) 

Sta. 1 


Suspension Line 
425 lb Cord 


Sta. 

AS (in. } 

Sw(in. ) 

27 

2.0 

5 . '5 

26 

2.286 

4.123 

25 

1.2 

5.022 

24 

0.0 

5.332 

23 

2.381 

5.947 

22 

0.0 

6 , 562 

21 

1.2 

6.871 

20 

3.81 

8.165 

19 

1.2 

9.458 

18 

0.0 

9.767 

17 

2.381 

10.382 

16 

0.0 

10.997 

15 

1.2 

11.307 

14 

3.81 

12.6 

13 

1.2 

13.893 

12 

0.0 

14.203 

11 

2.381 

14.818 

10 

0.0 

15.433 

9 

1.2 

15.742 

8 

3.81 

17.036 

7 

1.2 

18.329 

6 

0.0 

18.639 

5 

2.381 

19.253 

4 

0.0 

19.868 

3 

1.2 

20.178 

2 

3.06 

21.277 

LL 

2.0 >23.1 


Fig. 28 . Pilot Chute Structural Model 
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Fig. 29 • Canopy Pressure Distribution for the Pilot Parachute 
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skirt Distance from Skirt, in. 

Fig. 30. Sail Loading for the Pilot Chute 
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Fig . 31 . Radial Tape Loads for the Pilot Chute 


4.2.3 Pilot Parachute Film Analysis 

In order to verify the analysis used to predict the internal 
loads In Figures 30 and 31 , a comparison of predicted and 
measured canopy shapes Is made. The air-to-air and ground-to- 
air photographic coverage Is Inadequate for canopy measurements. 
Onboard coverage is used to compare canopy diameters as a function 
of the distance from the skirt. Figure 32 shows a comparison 
of predicted and measured canopy diameters for the pilot chute at 
a load of 3*600 pounds. The photograph used was taken just as 
the main parachute pack was lifted from the vehicle. Load traces 
indicate that the pilot chute had an axial load of 3*600 pounds 
at this time. The close correlation of the predicted and measured 
canopy shapes supports the internal load analysis of Section 3.1. 
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4.3 MAIN PARACHUTE 

4,3.1 Structural Models of the Main Parachute 


The construction of the main parachute requires an elaborate 
model to reasonably simulate the nonuniform geometry and stiff- 
ness. Such features as sail fullness and reinforcing bands are 
represented in t’.e model by horizontal elements with appropriate 
geometry and stiffness. Estimation of the width over which a 
heavy reinforcing member acts is problematical, but the overall 
solution is not overly sensitive to this modeling. Use of al- 
lowable strains rather than allowable load for stress analysis 
of the adjacent area gives a second check on the margins of 
safety. Slots are accounted for in the solution by traversing 
the slot segment with no change in meridional member load or 
slope. The canopy is reefed in two stages by reefing lines at- 
tached to the skirt band at midgore points. Main parachute con- 
struction is -hown in Figures 33 and 34 and the corresponding 
structural model is shown in Figure 35 . 

Three different models corresponding to the three major steps 
in the evolution of the final design are analyzed for the fully 
inflated condition. They differ only in the number and strength 
of reinforcing bands on the sails. 

Version 313 has additional tapes added to the trailing edges of 
sails, numbers 5 through 10. Version 519 has additional rip-stop 
bands at the trailing edges of Sails 5> 9 and 10. 

Solutions obtained for the three versions under identical loading 
conditions give an insight into the effect of the reinforcing 
bands on the internal load distribution. In Figure 36 , a 
comparison of the ratio of applied load to rated strength of the 
sail 3 for the three versions is given. It is seen that in ad- 
dition to providing additional strength locally, the bands re- 
duce the load in adjacent areas by reducing the overall diameter 
of the canopy. The critical areas which occur at the trailing 
edges of the sails with leading edge fullness give the highest 
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peaks. The trailing edges of Sails 7, 8, 9 and 10 are the 
critically loaded horizontal members. 

Sail 8 trailing edge has a load-to-strength ratio of 0.55 for 
the 507 version. If this member had an efficiency of 55 per- 
cent, a failure would be predicted at an axial load of 52,800 
pounds. The addition of reinforcing tapes to the parachute to 
form the 513 version reduced the load to strength ratio at the 
trailing edge of Sail 8 to 0.45. If this area remained critical 
in the parachute, the axial load capability of the 513 parachute 
would be increased by ten percent over the 507 version. How- 
ever, the trailing edge of Sail 9 is seen to be critical for the 
513 version with a load to strength ratio of 0.51. 

The load to strength ratio for Sail 9 in the ^19 version is re- 
duced to 0.24 by the addition of two 1200-pound tapes to the 
trailing edge. The failure analysis in Section 4.3.3 verifies 
the failure predictions of Figure 36 . 
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Fig. 33 Main Parachute Assembly 
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Fig* 35 • Main Parachute Structural Model 
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4.3-2 Main Parachute Failure Analysis 

Several versions of the R766I main parachute were tested during 
the evolution of the final design. For this analysis these 
versions are analyzed in three groups. The versions in each 
group are essentially the same structurally and the groups cor- 
respond to the three versions analyzed in Section 4.3.1 All ol 
these versions have 75 percent of Sr II 5 removed. 

The first generation (Group l), consisting of versions 501 through 
507 t had five ultimate strength tests: Drop Tests 29-3 A, 46-2, 

46-1A, 46-3 and 48-4. Group 2 consisting of versions 509 through 
513» had three tests at or near design ultimate load: Drop Tests 

46-4, 46-5 and 70-3. 

Group 3> which includes the present spacecraft configuration, 
consists of versions 515 through 519- This group has had three 
ultimate load tests in the full open condition: Drop Tests 

30-2, 30-3 and 82-4. 

Group 1 canopies were an early version with a minimum of sail 
edge reinforcements (See Figure 3^ )• Using the analysis of 
Section 4.3.1, together with appropriate joint efficiency and 
design factors (net efficiency = O.52), gives a predicted ulti- 
mate strength 31 ,100 lb, with failure starting at Sail 8. Of 
the five tests, three are considered nonrepresentative because 
of the following reasons: In Test 46-1 A, the parachute was 

overloaded to 33>900 lb in the reefed condition, which probably 
damaged critical joints. In Tests 46-3 and 48-4, Sail 8 failed 
by splitting along the gore centerline. Since the radial seam 
joint is normally weaker than the base material, this type of 
failure indicates that the sail cloth was damaged in this area. 
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Post test examination of the parachutes used in Tests 29-3A, 

46-1A and 46-3 revealed creases and small splits along the gore 
centerlines of several sails. This damage was attributed to 
improper packing procedures. In the remaining two tests of 
Group 1, canopy failures occurred at leads of 27,600 lb in 
Test 29-3A and 35 , 050 in Test 46-2. These values vary about 
13 percent from the predicted failure load of 31 ,100 lb. 

The first two tests of Group 2 canopies show good agreement 
with predicted strengtns. In Test 46-4 no major damage occurred 
at a load of 32,350 lb, and in Test 46-5 which reached a load 
of 35,060 lb, a gore split started at a Sail 7 radial seam. 

The predicted ultimate load is 33,^00 lb with failure at the 
trailing edge of Sail 9. The failure that occurred in Test 
70-3 at 22,390 lb cannot be explained by the stress analysis. 

Both Group I and Group 2 canopies had been subjected to higher 
loads than that which caused failure in Test 70-3- However, 
these higher loads were achieved by Increasing test dynamic pres- 
sure, whereas in Test 70-3 a heavier test vehicle was used to 
achieve the high loading. The effect of the heavy test vehicle 
is that peak loading occurs later in the opening process, there- 
by increasing the loading in the lower sails. In the post test 
evaluation of this test, it was concluded that the canopy was 
structurally inadequate, and as a result, heavy reinforcing 
rings were added to the trailing edges of Sails 5 and 9. The 
separation of the pilot chute during Test 70-3 was not considered 
to have been an important factor at that time. In two later tests, 
80-3 and 80-3R, pilot chute separation preceded failure by gore 
splitting at near limit load (reefed). Although the mechanics 
of the relationship cannot be determined with the information 
available, it is believed that these failures were related in 
Test 70-3. 
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Group 3 canopies have heavy "rip stop" bands on Sails 5 and 9 , 
which increase the load required to start gore splitting, and 
also prevent any splits which might result from deployment 
damage from propogating to canopy collapse. The predicted 
load for sail failure is approximately 40,000 lb for Group 3 
canopies, but the ultimate strength is limited to 33*700 lb 
by the strength of the suspension lines. This prediction is 
supported by the three ultimate strength tests of Group 3 para- 
chutes. In Test 30-2 no major damage occurred at 32,840 lbs. 

In Test 82-4 the three suspension line breaks which occurred at 
32,200 lb were attributed to deployment damage in post test 
evaluation. In Test 30-3 catastrophic failures of suspension 
lines and radial tapes occurred at 34,300 lbs. 

A summary of the main parachute failure analysis is presented 
in Table 2. 
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Table 2. Main Parachute Failure Analysis 
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29-3 A and 48-4 were two chute clu»* er teats. The unsysmetrlcal shapes developed due t> 
aerodynamic Interference may increase local loading by more than the 1.05 factor used : 
the analysis to account for unaymr itrloal loading. Lack of data from cluster testa at 
ultimate load precludes making a quantitative evaluation of this factor. 



4.3*3 Main Parachute Pressure Distributions and Film Analysis 

Pressure distribution data for reefed parachutes are not avail- 
able. In the stress analysis of the Apollo Earth Landing System, 
it was assumed that the pressure distribution at peak ri?er 
load was uniform. This study shows that pressure distributions 
can be approximated by comparing predicted and measured canopy 
profiles. 

A review of the Apollo Drop Test film coverage failed to produce 
adequate profile views of the reefed main parachutes at the time 
of peak loading. Two usable profiles were found, however for 
the parachute in steady state descent. 

A profile of the first stage reefed canopy was found in alr-to 
air coverage of Drop Test 46-1A. By using tne widths of the 
slot and a sail panel as reference lengths, dimensions of the 
canopy were calculated from the photographic image. In the 
first stage the canopy was observed to inflate from the vent 
to the Sail 5 slot only. Using the leading edge of Sail 4 as 
a starting point, the analysis method given in Section 3*1 was 
used to solve for the equilibrium canopy shape. A range of 
pressure distributions as shown in Figure 37 is used to de- 
termine the sensitivity of the canopy profile to variations 
in pressure. The first stage pressure distributions have zero 
pressure from skirt to the leading edge of Sail 4 (66 percent 
of the distance from skirt to vent), linearly increasing pres- 
sure to three locations above this point, and then constant 
pressure to the vent. Figure 38 shows the observed first 
stage profile compared to predicted profiles for the three 
first stage pressure distributions in Figure 37 . The constant 
pressure shape (distribution A) is seen to have a flatter pro- 
file than that observed in photographs. Pressure distribution C 
produces a profile that is in close agreement with the observed 
steady state profile. 
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NORTHROP 

A peak can be seen In the first stage profile at the vent. This 
appears to be due to the pull of the attached pilot chute. No 
provision has yet been made to the theoretical model to account 
for this affect. 

Since the canopy is not critical during steady state descent, 
the internal load predictions are not included here. For the 
internal load predictions of the main parachute during the open- 
ing process, see the stress-time study in Section 4.3.4. 

A profile of the second stage reefed main parachute was found 
on ground-to-air coverage of Drop Test 82-1R. The canopy was 
observed to inflate from the vent to the leading edge of Sail 8. 
Using this starting point, the analysis method given in Section 
3.1 was used to solve for equilibrium canopy shape. The second 
stage pressure distributions shown in Figure 37 are used to 
show the sensitivity of the canopy profile to variations in 
the pressure distribution. The differential pressures are zero 
below Sail 8 (33 percent of the distance from skirt to vent), 
increasing linearly to three locations above the first inflated 
sail, and then remain constant to the vent. Figure 39 shows 
the observed profile compared to predicted profiles for the 
three second stage pressure distributions in Figure 37 . Pres- 
sure distribution C produced best agreement with the measured 
profile but further refinement of the pressure distribution 
would improve this comparison. Since the steady state descent 
phase of the second stage is not critical, the Internal load 
predictions are not presented here. The internal loads at the 
critical times in the opening process are discussed in Section 
4.3.4. 

Pressure distribution for the fully Inflated main parachute is 
approximated by interpolating the wind tunnel test data taken 
from Reference 20 for a solid circular and a ringslot 
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parachute. Figure 40 shows the resulting distribution with 
peak pressure at the skirt. While the infinite mass condition 
of wind tunnel testing is considerably different from the drop 
test conditions being analyzed, this is the oest data available, 
and the analysis which follows shows reasonable correlation. 

It has been observed in drop test films that the parachute is 
not inflated to its maximum diameter at the time of peak riser 
load. (See stress-time analysis. Section 4.3.4). For this 
reason the disreefed canopy is analyzed at both the time of maxi- 
mum diameter and time of peak riser load. 

Figure 4l shows that the pressure distribution described above 
produces a good comparison of predicted and measured canopy shapes 
for the time of maximum diameter. The measured canopy shapes are 
from the onboard coverage of Drop Tests 30-2 and 82-4. The method 
of obtaining canopy dimensions from the photographic image is de- 
scribed in Section 4.1.4. Measurements of the canopy at the time 
of peak riser load, however, indicate a smaller canopy. By ex- 
tending the pressure distributions used In the analysis of reefed 
canopies, the second distribution in Figure 40 was found to 
produce a shape that agreed with drop test measurements as shown 
in Figure 42 . The pressure distribution indicates low pressure 
at the skirt, increasing linearly to a location above the skirt 
equal to 12 percent of the radial tape length, and then remaining 
constant to the vent. Internal load predictions for the canopy 
at the times of maximum diameter and peak riser load are presented 
in Section 4.3.4. 

Results of this study show that the differential pressure distri- 
butions on reefed and open canopies can be approximated from 
photographic records of canopy shape. 
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Fig. 37. Reefed Main Parachute Pressure Distributions 
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Comparison of Measured and Predicted Profiles for 
the First Stage Reefed Main Parachute 
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Fig. 39 . Comparison of Measured and Predicted Profiles for the Second Stage Reefed Main Parachute 
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Fig. 40. Fully Inflated Main Parachute Pressure Distributions 
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Fig. 4l . Comparison of Predicted and Measured Main Parachute 
Diameters at the Time of Maximum Canopy Diameter 
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Comparison of Predicted and Measured Main Parac 
Diameters at the Time of Maximum Riser Load 
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4.3.4 Stress-Time Study of the Apollo Main Parachute 

Previous analyses of the main parachute considered only three 
Instants of the Inflation process. First stage reefed, second 
stage reefed, and full open. It was assumed that maximum in- 
ternal loading occurred at the instant of maximum riser loading 
for each of the three stages and that all the sails in the 
stage were inflated at this time. The present analysis improves 
the margin of safety calculations by eliminating the above as- 
sumptions. 

In order to correlate analystlcal results with empirical evidence 
the canopy is analyzed at times corresponding to the instant 
each sail becomes inflated. The term "inflated" is defined as 
the time at which a sail ceases to flutter randomly. The 
empirical data needed includes the riser load and pressure dlstri 
bution over the canopy as each sail inflates. Photographic 
coverage provides the time at which each sail becomes inflated 
and the load-time trace gives the riser loading at each of 
these times. Since the parachutes were not instrumented to re- 
cord pressure distribution, the canopy shapes are used to infer 
the pressure distribution. 

Drop Test 82-2 was found to be the most acceptable test to pro- 
vide the needed empirical data. The test resulted in high canopy 
loadings for all three stages of inflation and had good cor- 
relation between film coverage and the load trace. A comparison 
of test loads with design ultimate loads is given below: 

Stage Design Ultimate Max. Test Load 

Stage 1 Reefed 29,700 lbs 20,375 lbs 

Stage 2 Reefed 32,000 lbs 32,710 lbs 

Stage 3 Disreefed 30,900 lbs 28,135 lbs 
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The test loads at the selected times are multiplied by the ratio 
of the design ultimate to the maximum test load for that stage. 
This calculation is valid because the internal loads vary linearly 
over short ranges of loading as was shown in Section 4.1.3. 

Table 3 gives the design loads used in this study. 

By using assumed pressure distributions and the design loads 
listed above, the internal loading at each state of inflation 
was calculated by computer program CANO. The computer program 
also calculated canopy shapes for each state of inflation. The 
analytical shapes were compared to measured shapes from photo- 
graphic coverage of Test 82-2. If the shapes did not agree for 
a particular state, the pressure distribution was revised and a 
new computation was made. This procedure was continued until 
the analytical canopy shapes at each state of inflation matched 
the measured shapes. Figure 43 shows the pressure distributions 
that produced good agreement of canopy shapes. 

The pressure distributions indicate that as the inflation process 
continues, the inflated portion of the canopy has low pressure 
at its leading edge, increasing to a uniform pressure over the 
upper portions of the canopy. The uninflated portion of the 
canopy is assumed to have zero differential pressure. For ex- 
ample, at the instant Sail 6 inflates, 49 percent of the canopy 
is un'nflated and the leading edge of Sail 6 has a differential 
pressure of 0.1, increasing linearly to a point 62.5 percent of 
the distance from the skirt to the vent. The pressure then re- 
mains uniform to the vent. At the time of maximum canopy dia- 
meter, the pressure is maximum at the skirt. Accurate profile 
views of the inflating canopy would improve the assumed pres- 
sure distributions shown in Figure 43. The best photographic 
coverage of Test 82-2 was provided by the ground-to-air cameras. 
Profile views of the canopy were limited to the early states of 
inflation. 
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The shape of the canopy after first stage disreef was based on 
measurements of maximum diameter, skirt diameter and suspension 
line angle. Some error can be exoected in the pressure distri- 
butions obtained by matching the analytical shape to only these 
three measured quantities. The actual canopy dimensions for 
each state of inflation were calculated by using the length of 
various canopy elements, such as slot width, sail '-'idth and 
reefing line diameter as reference lengths. 

Figure 44 shows the riser load-time trace for the design loads 
in Table 3. Figure 45 shows the critical sail loadings as a 
function of time. The critical sail loading for the entire 
inflation process occurs at Sail 1 just as Sail 4 inflates. 

The riser load at this time is 29,700 pounds as compared to a 
maximum riser load of 32,200 pounds occurring just as Sail 6 
inflates. Figures 46, 47 and 48 show the loading on each 
sail throughout the inflation process. 

The dashed lines in Figures 45, 46, 47 and 48 indicate assumed 
sail loadings for the portions of the inflation process that 
were not analyzed. These portions of the process are at times 
of decreasing riser load when the canopy is approaching a 
steady state condition. 

Results of this study show that the assumption that the maximum 
sail load-strength ratio for each stage occurs at the time of 
maximum riser load for that stage is valid. The assumption 
that maximum riser load for a stage occurs when all sails of 
that stage become inflated is valid for the first reefed stage 
and full open, but is not true for the second reefed stage. 

Peak riser load for the second reefed stage occurs just a3 
Sail 6 becomes Inflated, 1.08 seconds before complete inflation 
of this stage. 
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Table 3 . Load-Time History Used for the Stress-Time Study 


Inflation 

State 


Sail 2 
Inflated 

Sail 3 
Inflated 

Sail 4 
Inflated 



16,430 


17,970 


20,375 


Load 

Ratio 


29,700 


20,375 

29,700 


20,375 

29,700 


20,375 



M.C.D.R. 1 


Sail 6 
Inflated 


32,710 


32,200 


32,710 


Sail 7 
Inflated 

Sail 8 
Inflated 


26,780 


19,330 


32,200 


32,710 

32,200 


32,710 


35.27 


M.C.D.R. 2 

Sail 9 
Inflated 

Sail 10 
Inf] ated 

Sail 11 
Inflated 

Peak Load 


Sail 12 
Inflated 

Sail 12 
Expanded* 


18,140 


22,380 


26, 840 


28,135 


26,782 


15,370 


30, 900 


28,135 

30,900 


28,135 

30,900 


28,135 


30,900 



28,135 


30,900 


28,135 


* At time of maximum skirt diameter. 


NOTE: M.C.D.R. denotes "Main Canopy Disreef" 
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Fig.43. Canopy Pressure Distributions Used for Stress-Time Study 
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Time From Launch, sec 

Fig. 44. Load-Time Trace for the 83.5 ft. D 0 Ringsail Parachute 
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Time From Launch, sec 

Fig. 45 . Critical Sail Load vs. Time for the 83.5 ft D Ringsail Parachute 
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Time From Launch, sec 

Fig. 46 • Sail loading vs. Time for Saii3 1, 2, 3, 
4 of the 83.5 ft D 0 Ringsail Parachute 


NORTHROP 



108 


NVR-6432 


Time From Launch, sec 

Fig. ^7. Sail Load vs. Time for Sails 5, 6, 7, an 
8 of the 83.5 ft D Q Ringsail Parachute 
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Table 4 shows the maximum radial tape load-strength ratio for each 
of the states of Inflation. Peak radial load of 68.1 percent of 
the rated strength occurs Just as Sail No. 6 Inflates. Thus, the 
radial tape load-strength ratio Is maximum at the time of maximum 
riser load. 


Table 4. Load-Tlme History of Main Parachute Radial Tape 


Time from 

Number of 

Maximum Radial | 

Location of 

Launch 

Sails Inflated 

Tape Loading 

Maximum Load 

sec 

1 

Load 


i 

4 * 

Rated Strength 


23.94 

! 2 

} 

0.505 

; L.E., Sail 2 

24.09 

• ? 

0.553 

Middle, Sail 3 

24.30 

i » 

! 6 

0.631 

Middle, Sail 4 

27.90 

0.681 

Middle, Sail 6 

28.31 

1 

7 

0.564 

L.E., Sail 6 

28.98 

! 8 

0.416 

Middle, Sail 7 

34.74 


0.444 

Middle, Sail 8 

34.92 

10 

0.525 

' Middle, Sail 9 

35.11 

11 

0.634 ' 

' L.E., Sail 10 

35.36 

12 

0.635 

l 

L.E., Sail 11 

35.63 

12 

0.367 

T.E., Sail 12 


i . (Max. Diameter) 

l i L 


Margin of safety calculations for the critical entry condition are 
presented in Appendix B. 
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4.3.5 Optimum Weight Calculation for the Main Parachute 

Since weight is one of the most critical parameters In designing 
spacecraft recovery systems , an insight to the theoretical 
weight optimum of a parachute would be useful. A theoretically 
optimum parachute would have continuous members (100 percent 
joint efficiency) and materials that would be just strong enough 
at every point to withstand the design loads. 

Although this optimum weight is a value that can never be achieved 
in practice* it serves as a yardstick for judging the efficiency 
of a parachute design. Optimum weight could also serve as a 
sound basis for predicting the weight of a new design. Optimum 
weight factors for existing parachutes of given type and class 
could be computed and averaged to obtain a typical value. The 
product of this factor and the calculated optimum weight for a 
new design of the same type and class would give a predicted 
weight. 

In this study* two optimum weight factors are calculated. The 
first factor is the ratio of spacecraft main parachute weight 
to optimum parachute weight including typical joint efficiencies. 

The second factor is the ratio of present main parachute weight 
to optimum parachute weight assuming 100 percent Joint efficiencies. 
This factor gives an absolute optimum weight by assuming all the 
members to be continuous as well as being of minimum strength. 
Redundant members (the two rip-stop bands and one of the first 
stage reefing lines) are omitted from the calculation* as are 
nonstructural items such as reefing line cutters and pockets. 

The strength-to-weight ratio for a fabric material is found 
by dividing its rated strength by its density in consistent 
units. The strength-to-weight ratio has units of length and 
can be visualized as the length that could just support its 
own weight. Hie values used in this analysis are 1.11 x lO^in. 
for tapes* 1.56 x 10° in. for cord, and 0.8l x 10° in. for cloth. 
These values are average minimums for fabric materials in the 
"Parachute Handbook". 
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Using the results of the parachute opening stress-time study 
(Section 4.3.4), an evaluation of optimum canopy weight was made 
by reducing the rated streng of each canopy member to giv' 
a zero margin of safety. Ef xciency factors such as joint ef- 
ficiency, abrasion, and unsymmetrical loading were included to 
give a realistic weight comparison with the present spacecraft 
parachute (thermal factors were neglected to simplify the 
analysis). After the rated strengths had been reduced, ad- 
ditional computer runs were made at each of the critical loading 
conditions using averaged load-strain curves. The equilibrium 
shape and internal loading was found for the pressure distri- 
butions shown in Figure M 3 . Based on the results of these 
computer runs, the rated strengths were again revised to give 
zero margins of safety for all elements and further computer 
runs were made. This procedure was continued until the computed 
internal loads resulted in zero margins of safety for all canopy 
elements. Average streng th-to-weight ratios for tape, cloth and 
cord were then used to compute the optimum canopy weight of the 
parachute. By using the same average strength-to-weight ratios, 
the weight of the spacecraft parachute was also found. The ratio 
of vne spacecraft main parachute weight to the optimum parachute 
weight with typical joint efficiencies is 1.53* To obtain the 
second of the optimum weight factors, the rated strength of each 
; lement In the parachute was reduced to give a zero margin of 
safety for 100 percent joint efficiency. The weight of the 
parachute was then found by using the same average strength-to- 
weight ratios as were used above. Table 5 contains the weight 
of each parachute component for the spacecraft design, for the 
optimum design with typical joint efficiencies and for the 
optimum design with 100 percent joint efficiencies. The ratio 
of spacecraft main parachute weight to optimum parachute weight 
with 300 ^ercent joint efficiencies is 2.17. 
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Table 5 . Comparison of Apollo Main Parachute Weight 
with Theoretically Optimum Weights. 


Item 

Spacecraft 

Design 

Optimum 

1* 

Optimum 

2** 

Suspension 

Lines 

42.10 lbs 

41.19 lbs 

30.89 lbs 

Sail 1 

1.15 

1.08 

0.54 

Sail 2 

2.53 

2.22 

1.11 

Sail 3 

3.89 

2.43 

1.21 

Sail 4 

5.38 

2.49 

1.25 

Sail 5 

3.29 

0.65 

0.34 

Sail 6 

4.14 

1.92 

0.97 

Sail 7 

4.77 

2.24 

1.11 

Sail 8 

5.13 

2.59 

1.30 

Sail 9 

9.72 

2.89 

1.44 

Sail 10 

5.30 

1.50 

0.75 

Sail 11 

3.75 

1.38 

0.69 

Sail 12 

6.12 

0.73 

O.36 

Skirt Band 

6.51 

3.00 

2.73 

Vent Band 

0.67 

0.56 

0.51 

Vent Lines 

0.62 

0.26 

0.16 

First Stage 
Reefing Line 

0.52 

0.15 

0.14 

Second Stage 
Reefing Line 

0.89 

0.49 

0.33 

Radial Tapes 

22.00 

17.66 

14.30 

Total 

130.48 lbs 

85.43 lbs 

60.13 lbs 


* Optimum 1 Includes typical Joint efflclences. 
** Optimum 2 Includes 100$ Joint efficiencies. 
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SECTION 5.0 

PILOT CHUTE RISER DYNAMICS 

A series of riser dynamics tests were performed as part of an 
effort to determine the cause of two failures that had occurred 
in Apollo drop tests conducted in August of 1967 (Tests 80-3 
and 80-3R). In both drop test cases, the pilot chute riser 
failed at time of main parachute canopy stretch. Eoth failures 
were, therefore, due to impact loadings rather than parachute 
opening loads. The predicted loads, using analysis methods 
employed at that time, were well below the strength of the riser 
assembly. In addition, the failures occurred in free lengths, 
rather than in the joints which are weaker than the parent 
material. In order to study the mode of failure, the fabric 
riser dynamics tests were conducted in the Northrop Ventura 
laboratories. The resulting data are the subject of this 
analysis . 

Longitudinal impact loading was applied at a velocity of 
300 ft/sec to samples of pilot chute riser material (nylon 
webbing) and structure elements. Test data consisted of high 
speed film sequences and strip chart force gage records. It is 
the purpose of this study to explain these data and to prepare 
recommendations for additional riser dynamics test and analysis. 

The approach employed in the endeavor to understand the phenomena 
occurring during the impact loading of the specimens is to first 
perform a simplified analysis of the observed wave motion. The 
method of characteristics is applied to an approximate mathematical 
model obtained by the application of eleven assumptions regarding 
the system boundary conditions, material properties and general 
behavior. The most significant of these approximations are the 
assumptions that the nylon webbing used in the specimens behaves 
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in a linear-elastic manner and that the elastic modulus is 
independent of loading rate. The application of this method of 
analysis contributes significantly to the understanding of the 
performance of the test specimens but is unsuitable for the 
prediction of webbing stresses; l.e., the assumptions used are 
excessively restrictive. This work is described in Section 5.4. 

When the complexity of the problem is increased by the necessary 
annulment of several of the assumptions listed in Section 5.4, 
the use of a high speed computer becomes desirable. Section 5.5 
is a description of a computerized application of the finite 
difference method of solution to the system of partial differential 
equations, boundary conditions and Initial conditions appropriate 
to more accurate representation of the real system. 

The arrangement of test apparatus used in these tests is shown 
in Section 5.1 and descriptions of the test specimens and 
relevant material properties are presented in Section 5.2. Forcq 
gage data are presented in Section 5.3. Seven film sequences 
have been studied quantitatively and plots of position versus 
time have been constructed for each specimen gage mark. A pre- 
sentation and interpretation of these results is given in 
Section 5.6. 

5.1 TEST EQUIPMENT AND PROCEDURE 

The equipment used in the riser dynamic tests is shown in 
Figure 49 . 

The specimens were suspended under light tension between the 
loop release mechanism and the resistance element type force 
gage. The force gage was attached to an extremely rigid Mad- 
man beam. Gage marks were painted on the trunk of the spe imen 
at two-foot intervals. 
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Pig. 49 . Pilot Chute Riser Dynamic Test Apparatus 



The impact wheel was brought up to a rotational speed that gave 
a rim speed equal to approximately 300 ft/sec. The impact loop 
was then released and allowed to be engaged by the wheel rim 
hcok. The resulting motion of the system was photographed 
against a background grid using high speed motion picture cameras. 
The force gage output was preserved by means of a strip chart 
recorder. 

The (nominally) ten-inch diameter impact loop was fabricated 
from 6000-pound test tubular nylon cord using a "Chinese finger" 
Joint with sufficient overlap to provide four active load carrying 
cords. There was no quality control on the fabrication of this 
part of the test structure. 

The high speed films are of variable quality. Only the tests 
for which the high speed films are readable are considered in 
this report. These tests and the respective impact velocities 
are listed in Table 6 . 

S2 DESCRIPTION OF SPECIMENS AND MATERIAL PROPERTIES 

The construction and materials of the nylon webbing specimens 
are shown in Figures 50 and 51. All material of a given speci- 
fication was obtained from the same lot. Notice that one of 
the webbings in the trunk of specimen type D4 was deliberately 
made approximately 3 percent longer than the other. The purpose 
of this was to determine the effect of mismatching on riser 
strength. 

The unstrained lineal density and the static load-strain curve 
for each type of webbing used in the test specimens are given 
in Figures 52 and 53. The straight line approximation used in 
the approximate linear-elastic analysis are superimposed over 
the nonlinear static curves. 
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Drid’’ e Trunk 

in. ID \ MIL-W-4088 1 in. x MIL-W-5625 3/4 in.-r"^ \ ,-3/8 in. ID 
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Fig. 50 . Pilot Chute Riser Test Specimens D3> D4 and d6 
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Pig. 51 . Pilot Chute Riser Material Test Specimen 



Perce, lb. 



Fig. 52 . Static Load - Strain Curve for a Sample 
of MIL-W-5625 3/4" Nylon Webbing 
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Load, lb 




8000 


6000 


4000 


2000 



Strain 

Pig. 53 . Average Static Load - Strain Curve for 
MIL-W-4088 1" Nylon Webbing 
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MMINROP 

In order to use the nonlinear curves In a computer calculation. 

It is necessary to either store the curve in digital form or in 
the form of an approximating function. An approximating function 
in the form of a least squares best fit fifth order polynomial 
has been obtained for the MIL-W-5625 3/4 inch nylon webbing used 
in all four types of test specimens. The expression, plotted 
in Figure 52 is 

r ■» (1.827041 x 10 4 ) e - (3.550943 x 10 5 ) e 2 

+ (3.680853 x 10 6 ) e 3 - (1.164287 x 10 7 ) 

+ (8.459224 x 10 6 ) e 5 

Because of similarity of load-strain curves, this expression is 
used by scaling for the MIL-W-4088 webbing in Section 5-5- 

References 23 and 2 5 show that for the case of constant average 
strain rate, the load-strain curve of nylon yarn depends upon 
the rate of loading. Therefore, a priori selection of the 
static curves for use as dynamic curves in the case of nylon 
in the form of webbing is not justified. However, in Section 5.5 
the static curves are used as an approximation. 

It is conceivable that the correct dynamic equation of state for 
nylon webbing is a differential equation involving stress, strain, 
strain rat^ and perhaps even strain history^? These tests were 
not comprehensive enough to derive such an equation of state but 
in Section 5*6 an estimate of the load-strain path actually 
followed by the webbing in each test is given. 
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5.3 FORCE GAGE DATA 


Deadman force as a function of time Is shown for each test In 
Figures 54 - 60. The experimental data did not yield a force 
gage time reference with respect to other events of the test. 
Therefore, the zero time points on the horizontal axes are at 
arbitrary locations. The peak deadman force for each test is 
(riven in Table 6 . 

Inspection of the force gage curves shows a significant 
difference between the general shapes for the cases of single 
and double webbing trunks. For the single webbing case (Fig- 
ure 59 ) the initial part of the curve is sloped, with small 
oscillations, whereas for the double webbing case (all other 
50 ts) the initial part of the curve displays large oscillations. 
In both cases the initial oscillations tend to decay, and a 
steady positive slope occurs prior to failure. 

Table 6 shows that, on the average, the mismatched specimens 
suffered a 9 percent decrease in peak load relative to the 
correctely fabricated specimens. The average peak load 
attained by the simulated risers (D3 and d 4) is 15 percent 
greater than the peak load attained by the riser sample (D 6 ). 

5.4 THE LINEAR-ELASTIC APPROXIMATION 

5.4.1 Assumptions 

As a first step in developing an understanding of the performance 
of the specimens in the riser dynamics tests, the linear-elastic 
case for longitudinal waves is considered. Eleven assumptions 
are involvedin the construction of the mathematical model. 

They are as follows: 
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Force, lb 



0 2 4 6 8 10 12 
Time, milliseconds 

Fig. 5 4 . Deadman Force Gage Record for 
Pilot Chute Riser Test D3-1 
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Force, lb 



0 2 4 6 8 10 12 
Time, milliseconds 

Fig. 5 5. Deadman Force Gage Record for 
Pilot Chute Riser Test D3-2 
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Fig. 5 6 • De adman Force Gage Record for 
Pilot Chute Riser Test D3-3 
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Fig. 57 . Deadman Force Gage Record for 
Pilot Chute Riser Test D4-5 
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Force, lb 




Fig. 58 • Deadman Force Gage Record for 
Pilot Chute Riser Test D4-6 
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Fig. 59 . Deadman Force Gage Record for 
Pilot Chute Riser Test D5-1 
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Fig. 60 . Deadman Force Gage Record for 
pilot Chute Riser Test D6-1 
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a) The load-strain curve of the specimen materials under 
dynamic loading is independent of rate of loading and 
strain history and is related to the static load-strain 
curve in the secant fashion illustrated in Reference 27 
and applied to the materials of the subject tests as 
shown in Figures 32 and 53. 

b) The end of the specimen moves at the same velocity 
as the rim of the impact wheel once the hook is en- 
gaged. The impact loop is assumed to possess zero 
stiffness until it is extended into a straight line 
and is assumed to be perfectly rigid thereafter. 

c) The impact wheel rim speed remains constant at its 
initial value V 0 . (in reality, the wheel rim speed 
decreases due to the reaction of the specimen.) 

d) No material is wrapped on to the impact wheel. The 
relaxation of this assumption makes the problem non- 
linear even though assumption a) is maintained. 

e) Specimen displacements are small and lineal. 

f ) The wave propagation speed is very large relative 
to the velocity of the mass particles. Therefore, 
the velocity of the waves relative to the laboratory 
can be taken to be the same as the velocity in a 
reference system fixed in the specimen. 

g ) Variations in lineal density and stiffness occurring 
over distances of four inches or less are modeled as 
a single step. This assumption has the effect of 
eliminating the details of the construction of the 
specimens at the ends of the overlay reinforcements. 

h ) The stress, strain, displacement and velocity are 
uniformly distributed over the cross section of 
the specimen, even near discontinuities. 
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i) The specimen is perfectly fixed at the deadman end. 

J) There is no interaction between transverse and 
longitudinal waves. 

k) The kinetic energy due to cross axis strain is 
negligible. 

5.4.2 Method Of Analysis 

The behavior of a system to which the preceding assumptions 
apply can be predicted by means of the application of four 
basic results. 

a) The longitudinal wave propagation velocity in 
section j of the specimen is given by: 



where kj is the slope of the linear load-strain 
curve for the specimen material in section 
j and Pj is the respective lineal density. 
(Reference 24) 

b) Immediately following the initiation of the step 
velocity impact at t = 0, a strain wave propagates 
away from the impact end with constant velocity 

•> -W 

The strain in front of 
strain behind the wave 

e = Yo 
° C 1 
(Reference 24) 


the wave is zero and the 
is given by 
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c) When a strain wave impinges on a discontinuity in the 
specimen, it is, in general, attcnm-.nad or amplified 
and a reflection wave is generated. The reflection 
coefficient of the discontinuity between section j and 
section j + 1 is given by: 



and the transmission coefficient is given by 

T, = — (l - R-t ) (64) 

« c b J 

where a = j and b = j + 1 if the wave approaches 
from section j and a =■ j + 1 and b = J if the wave 
approaches from section j + 1. For a fixed boundary 
and the impacted end of the specimen, Rj =1. 

(Reference 28 ) 

d) The strain at any point in the specimen at time t is 
equal to the algebraic sum of the amplitudes of all 
strain waves that have passed that point in the time 
interval 0, t (Superposition principle). 

Using the above four results as rules of operation, the response 
of a specimen to the velocity impact can be constructed graphically. 
This has been done for the D5 test and the construction, which 
is an application of the method of characteristics, is displayed 
in Figure 6l. 
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Fig. 6l . Wove Diagram, Linear-Elastic 
Approximation, Test Type D5 
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In this wave diagram, the horizontal axis Is dimensionless 
distance from the Impact end, £ . The vertical axis Is dimen- 
sionless time T = t/tp where tp is the time required for a wave 
to propagate from one end of the specimen to the other. The 
slanting lines indicate the path of the wave fronts in space- 
time and the numbers associated with each of those lines are 
the amplitude of the respective wave front expressed In terms 
the dimensionless strain, 


where e is the local strain . 

The procedure used is to lay out the initial wave front locus 
(dimensionless slope = c-^ tp/L, dimensionless amplitude = 1) 
extending from the origin to the first discontinuity at which 
point the reflected and transmitted wave amplitudes are calcu- 
lated using Equations 63 and 64 The transmitted wave moves up- 
ward and to the right (in space-time) with slope corresponding 
to the wave propagation velocity in the .second section (dimen- 
sionless slope = Cp tp/L). The reflected wave moves upward 
and to the left back into the region corresponding to the first 
section. These two new waves will then enter into further re- 
flections at discontinuities and boundaries. As the construction 
proceeds, the number of waves multiplies at an increasingly rapid 
rate. In order to make graphical calculations feasible, waves 
with amplitude less than .03 are terminated and waves approaching 
each other within At = 0.05 are superposed. The strain level 
at any value of x/L and dimensionless time t is obtained by 
summing the amplitudes of all waves intersecting a straight 
line between the points (x/L,0) and (x/L, t). 
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5 . 4 • 3 D escription Of The Approximate Model 

The graphical wave tracing technique Is applied to Test D5-1, 
the specimen of which Is described In Figure 6l . The parameters 
of each section of the specimen are listed in Table 7 . If more 
than one discontinuity appears in any 4-ir.ch length of specimen, 
they are lumped together in^o one discontinuity in this table. 

5.4.4 Results 

The longest uniform section of the specimen, which is referred 
to as the trunk of the specimen, consists of webbings with a 
combined rated strength equal to or less than half of the rated 
strength of any other section in the specimen. The strain distri- 
bution in the trunk of specimen D5 near the experimentally ob- 
served time of failure is shown in Figure 62 . A significant 
observation is that at this instant in time, the maximum strain 
in the trunk occurs at a station away from either end of the 
trunk. This sample situation, therefore, implies the possibility 
of failure in such structures in free lengths of the structure 
if the stress concentrations at discontinuities are sufficiently 
mild. 

In Figure 63 » strain is shown as a function of time for both 
ends of the trunk and for the end of the webbing adjacent to 
the deadman. The sweep of the strain wave along the specimen 
is manifested in the sequence of the initial response at each 
station. Also shown in Figure 63 are the dimensionless ultimate 
strains for the material of the trunk of the specimen of Test 
D5-1 and one other D5 type test. There are twc values because 
the impact velocities were not the same in the two tests per- 
formed on D5 specimens. It is interesting to note that the 
ultimate strain level is penetrated almost simultaneously at 
opposite ends of the trunk. The test films (Section 5.6) reveal 
simultaneous failures at the same points at approximately the 
same time as indicated by this approximate analysis. 
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Table 7 . Specimen Parameters For Approximate 
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Fig. 62. Strain Distribution Along the Length of the Trunk 
of Specimen Type D5. Approximate Linear-Elastic 
Analysis, T =* 2.8. 
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Fig. b3. Strain as a Function of Time at Three Stations, 
Approximate Linear-Elastic Analysis 



A general observation of importance is that at all of the three 
stations shown, the webbing is subjected to several cycles of 
loading and unloading before failure occurs. 

At the fixed end (deadman) of the specimen, the strain levels 
are generally lower than in the trunk because of the multiple 
layers of webbing used for reinforcement. In the tests, the 
load at the deadman was measured by a force gage and a comparison 
with the approximate analytical result is appropriate . Figure 
Ck gives the comparison for the case of Test D5-1. It is clear 
from this illustration that among the list of assumptions in 
Section 5.4.1 there is some combination that contributes to 
the introduction of significant errors. 

An observation of interest is that the calculated oscillations 
in the strain in the specimen at the deadman end are of approxi- 
mately the same period as oscillations found in the force traces. 
It can be seen from the wave plotting diagram that these oscil- 
lations are, at least In part, due to repeated reflections be- 
tween the opposite ends of the overlay reinforcement at the 
deadman end of the specimen. It is also apparent that the 
deadman force-time history has a first order dependence on the 
nature of all discontinuities, even those at stations remote 
from the deadman end. 

If wave propagation effects are neglected, that is, the tensile 
force in the specimen Is assumed the same at all stations, then 
the deadman force as a function of time would be represented by 
the straight line in Figure 64 • It can be seen that this repre- 
sentation is not adequate and that wave propagation effects must 
be Included for accurate structural analyses of pilot parachute 
risers. 
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Deadman, Force, 1000 lb 



Time, milliseconds 


Fig. 64 . Deadman Force as a Function of Time. 
Specimen Type D5. 
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5.^.5 Effect Of Finer Resolution Of Discontinuities 


In order to determine the effect of including the details of 
the construction of the outboard end of the deadman end 
reinforcement in Specimen Type D5* another wave diagram was 
constructed. The parameters used in this fine lumping case 
are listed in Table 8 . The results of this calculation are 
compared with the early time results (t <1.2) of the previous 
calculation in Figure 65. Examination of this figure shows 
that the effect of the finer resolution is to decompose each 
step into a sequence of smaller steps and to increase the 
magnitude of the initial jump in strain. The magnitude of 
these effects is much less than the magnitude of discrepancies 
in the previous results. Therefore, the lumping of discon- 
tinuities Is not considered a significant source of error. 

5.4.6 Observations Regarding The Approximate Analysis 

a) The discontinuities in the specimen affect the 
strain history at all stations in che specimen. 

b) All stations in the specimen are subjected to 
cycles of loading and unloading. Therefore, 
hysteresis effects must be considered. 

c) The stress at some station away from the ends of 
a uniform section of webbing may become greater 
than the stress at either end. 

d) Multiple failures are to be expected in destructive 
tests. 

e) Wave phenomena must be considered in the analysis 
of the riser dynamics tests. The uniform force 
assumption is not suitable. 

f) This approximate model is not an adequate 
representation of the real system. 
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Table 8 . Parameters for Fine Resolution of 

Discontinuities in Specimen Type D5 


Section 

Left 

Station 

(in.) 

Right 
Station 
(in .) 

Lineal 

Density 

(sl/ft) 

Stiffness 

(lb) 

Wave 

Velocity 
(ft/sec ) 

1 

0.00 

95.00 

.00115 

24,800 

4640 

2 

95.00 

263.50 

.000574 

12,400 

4640 

3 

263.50 

264.50 

.00115 

24,800 

4640 

4 

264.50 

265.50 

. 00172 

37,200 

4640 

5 

265.50 

276.00 

. 00230 

49,600 

4640 
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Pig. 65 . Strain at Deadman End of Specimen for the Cases of 
Coarse and Pine Resolution of Discontinuities 




5.5 The Finite Difference Method 


Section 5 . ^ contains a list of assumptions which, when applied 
to the riser dynamics test analysis, leads to unsatisfactory results. 
If these assumptions are revoked, the complexity of the resulting 
analysis increases to a level at which the use of high speed 
digital computers becomes a necessity. 

Digital computers may be applied to this problem in three ways. 

(a) Integration of the equations of motion of an 
approximate system consisting of a string of mass 
elements connected by suitable force generating 
elements, which would be linear springs in the 
linear-elastic case. 

(b) Application of finite difference techniques to the 
direct integration of the system of equations of 
motion, continuity equations, equations of state 
and boundary and Initial conditions. 

. 29-31 

(c) Application of the method of characteristics. 

Methods (a ' and(b) are equivalent and possess the disadvantage 
that spurious oscillations are introduced into the wave profiles. 
Reference 33 presents a modification of method b) in which the 
concept of artificial viscosity is introduced to suppress these 
oscillations. The availability of a computer program (AFTON) 
which implements this modification of method b) led to the 
selection of this method as the analytical procedure of analysis 
to be used in this study. 

An area of progress in this effort has been in the study and 
modification of the AFTON program. This program, which is 
described in Reference 34 , computes longitudinal wave motion 
in one -dimensional continua by solving the finite difference 
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forms the appropriate hydrodynamic equations. Several 
required program modifications have been made. 

The program was updated so that it would run on the Northrop 
IBM 360/65 computer. The logic was corrected so that the program 
would properly calculate the motions of the type observed in the 
laboratory tests. Also, extensive modification was completed 
that would allow simulation of the elimination of riser material 
as it passed on to the impact wheel rim. Within the duration 
of a typical laboratory test, a large protion of the riser 
sample passed onto the rim. Therefore waves must be made to 
reflect from the rim of the flywheel rather than from the loop 
end of the riser. The program logic has been changed so that 
it can accommodate this unusual situation. 

A shortcoming of the finite difference method is that, in order 
to avoid the fluctuations behind a wave front as observed in the 
lineal spring -mass model, an "artificial viscosity" is introduced 
into the computation. This artificial viscosity has the (un- 
desirable) effect of spreading an infinitesimally thin wave 
front (as occurs in the elastic case) into a computed wave front 
of finite thickness. In other words, after the continuum has 
been discretized into zones, the artificial viscosity causes 
wave fronts to be smeared over the length of several zones. 

This is undesirable because there were processes, sometimes im- 
portant, which took place over lengths much shorter than program 
zone lengths. Thus, the situation is as follows. In order to 
reduce or eliminate the fluctuations after a wave front has 
passed, the artificial viscosity must be increased. But in- 
creasing the artificial viscosity smears the wave front over 
a greater number of zones, causing the loss of fidelity of 
processes over short lengths. This fidelity can be improved 
by increasing the number of zones in a given length, thereby 
reducing the actual length over which the wave front is smeared. 
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But such an approach has the disadvantage that the computation 
time increases almost proportionately with the square of the 
number of zones. So an attempt has been made to assess the 
trade-offs between residual fluctuations, wave front smearing 
and computation time; and to choose values of l) the number 
of zones and 2) the Richtmyer-von Neumann^ artificial visco- 
sity constant that give reasonable results in terms of problem 
requirements. 

An adequate number of combinations of zone length and artificial 
viscosity were tried to allow their respective influences to be 
estimated for the range studied. Stated otherwise, given 
limits on load overshoot and wave front smearing, satisfactory 
values of zone length and artificial viscosity can be chosen, 
and computer time can be estimated. Unfortunately, the shortest 
lengths (1 in.), which occur in the tapered joints, and which 
are believed to be important, cannot be modeled in the AFTON 
program because they are too short. Therefore, it will probably 
have to suffice to include the effects of the next shortest 
pieces, which are 1 ft in length. 

Figure 66 is included to Illustrate a typical calculation done 
by the nFTON program. The trace represents the force wave profile 
at one instant in tiiae. This particular calculation is for the 
combination of 398 elements in a 23-foot uniform elastic sample 
and an artificial viscosity constant of 2.56. It shows both the 
overshoot and wave smearing described above. 

The AFTON program was used to calculate the force gage response 
of Tests D3-2 and D5-1 for the case in which several of the 
linear-elastic approximate case assumptions listed in Section 
5.1 are revoked. The changes with respect to the linear-elastic 
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Pig. 66 • Wave Profile Determined by the AFTON 

Program for the Case of a Uniform Specimen, 
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approximate case that were made are listed below using the 
respective Section 5.4.1 identification letters. 

a) The material is assumed to follow the non- 
linear static load-strain curve in an elastic 
manner 

b) No change 

c ) The impact wheel rim speed is allowed to decay 
due to the reaction force of the specimen 

d ) The bridle is wrapped on to the wheel rim 

e) Large displacements are allowed 

f) No restriction on magnitude of wave propagation 
velocity relative to particle velocity 

g ) No change 
h-k) No change 

Figures 67 and 68 are comparisons of the calculated force 
traces with the experimental curves. Since the zero time 
reference of the force gage data is unknown, the curves are 
Indexed with respect to the first significant peak. Due to 
the presence of a short concave downwards portion near the 
origin of the load-strain curves, the computed force traces 
display a "precursor wave" effect prior to the arrival of the 
sharp edged wave associated with the concave upwards portion 
of the load-strain curves. Some of the experimental curves in 
Section 5.3 possess this feature but it is not as pronounced 
as in the calculated curve, perhaps due to viscoelastic 
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effects. In the case of test D3-2, following the precursor, 
the experimental and calculated curves have roughly the same 
general shape but It is clear that among the assumptions 
remaining, there exists some combination that is not justified. 

In the case of Test Dp-1, the initial oscillations following 
the precursor effect are more faithfully represented than in 
the case of the linear-elastic approximation (Fig. t>4 ) and the 
later portion of the force trace Is represented approximately 
as well if the time scale for the calculated curve is con- 
tracted by a factor within the range corresponding to the 
Intertest variations In wave velocity reported in Section 5.6. 
However, the results for Test D5-1 confirm the need for an 
Improved mathematical model. 
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5.6 ANALYSIS OF FILMS 


Seven film sequences have been analyzed. Section 5.6.1 Is 
an in d'pth review of one film sequence and Section 5.6.2 
Is a comparison of wave velocities and webbing strains as 
measured from all seven film sequences. 

Reference 35 shows that for some materials with nonlinear 
load-strain diagrams the strain wave front will have a non- 
zero thickness and will change in shape as it progresses along 
the specimen, that is, each strain level may propagate with 
a different velocity. If this type of motion existed in the 
riser dynamics tests, it was not detected in the films due to 
the level of resolution imposed by the film framing rate. 
Therefore, the wave velocities obtained from reduction of the 
film data are referred to in the following as "apparent wave 
velocities" because they describe a gross effect rather than 
the fine details of the motion. 

5.6.1 Analysis Of Test No. D5-1 

Longitudinal Waves, Particle Velocities, Strain 

Figure 69 is a plot in position-time space of the events 
occurring prior to failure In Test No. D5-1. The horizontal 
axis is distance from the deadman support and the vertical 
axis is time from the arbitrary beginning of the film sequence. 
The left boundary represents the point at which the bridle 
of the specimen is feeding on to the rim of the wheel. The 
right boundary represents the point at which the webbing is 
attached to the force gage link. 
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Pig. 69 . Space-Time Diagram, Test D 5 






The solid lines show the paths of the gage marks. Gaps in- 
dicate loss of data due to obscuration of gage marks by twisting 
of the specimen or poor contrast with the background. Dis- 
continuities in the slope of these lines indicate th ' passage 
of a strain wave. The locus of the leading sweep of the strain 
wave along the length of the specimen has been drawn in as the 
line AB which is judged to be the straight line which passes 
most closely through all of the initial major discontinuities 
in the gage mark curves. The slope of line AB indicates an 
apparent wave propagation speed, of 3790 ft/sec. 

The Intersection of this straight line with the left boundary 
of the plot gives an estimate for the time at which the wheel 
hook engaged the loop at the end of tne specimen. In this 
case t c - 0.9 milliseconds. A similar straight line BC drawn 
through the second major slope discontinuities in the loci of 
the gage marks reveals a reflected wave traveling from the 
rigid support back towards the wheel rim. The slope of this 
line indicates a wave propagation speed, (c^)^, of -4630 ft/ 
sec relative to the laboratory. 

Prior to the arrival of the wave AB, the gage mark loci are 
vertical, indicating that the webbing is stationary until 
the wave arrives. The action of the passage of the wave AB 
is to introduce a nearly constant velocity as indicated by 
the nearly constant slope in the gage mark loci between lines 
AB and BC. The velocity of the webbing relative to the 
laboratory indicated by the best fit slopes between lines AB 
and BC are given in Table 9. 

The velocities of gage marks 1, 2, 3» 4, 9 differ from the 
wheel impact velocity of -275 ft/sec by less than 7 percent. 

At gage number 7, the velocity is considerably different 

from the impact velocity, probably due in part to the comparable 
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ie 9. 


Gage Mark Velocii 
First Wave, Test 


Gage Mark Number 


1 




3 

4 


c 

s 


6 

7 

Average 


U\ 


After Passage of 

-1 


Velocity 

(ft/sec) 

-294 

-294 

-294 

-294 

-272 

-200 

-275 
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magnitudes of the time between the passage of the two waves 
and the tim^ between exposures in the test films which gives 
an insufficient sampling for the accurate determination of 
velocity at gage number 7. 

The average strains between gage marks between the passage of 
waves AB and BC can be computed from the distance between gage 
marks before and after th^ passage of wave AB and are tabulated 
in Table 10. The strain in the specimen after the passage of 
the second wave and the strain at the time of the first failure 
are also shown in Table 10. 

For the webbing material on which the gage marks are placed, 
the ultimate static strain xs given as 0.23 (Figure 52). The 
maximum strain discovered between the waves AB and BC is 57 per- 
cent of this value, after the wave BC it is 32 percent and at 
the time of failure it is 90 percent. 

It is of interest to compare the velocity of the waves AB 

and BC relative to a coordinate system fixed to the specimen 

(the "Lagrangian velocity") with the velocity predicted by 

the linear elastic analysis. Jsing the results of Fenstermaker 
•32 

and Smith J , the Lagrangian velocity c of a longitudinal wave 
Is given by 

c = 1 + e ( c ," v ) 

where e is the strain in the material into which the wave is 
traveling, c is the wave velocity with respect to the laboratory 
and V is the velocity of the material into which the wave Is 
propagating. Using as V the average gage mark velocity after 
the first passage of the wave and as e the respective average 
of the strains between the gage marks, *• ^parent Lagrangian 

velocities are: 

fc >ab = rhj ( 379 ° + 0 ) - ~ lt/3ec 

(0> BC = I v( 09 06 (' 463 ° + 275 ) ’ - 3990 ft/3eC 
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Table 10. Webbing Strain, Test D5-1 


Between Gage 
Marks . . . 

Strain 

Between First 
Two Waves 

After Second 
Wave 

At Time of 
First Failure 

1 and 2 

.0924 

.109 

.160 

2 and 3 

.1304 

.189 

.207 

3 and 4 

.0598 

.128 

.179 

4 and 5 

.0796 

.150 

.186 

5 and 6 


• 064 

0 ± 

1 
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The theoretical wave velocity calculated from the "secant" 
approximation used in the linear-elastic analysis gives a 
Lagranglan velocity equal to 4650 ft/sec for both waves. The 
first wave AB propagates with a velocity that is 82 percent of 
the linear-elastic case and the reflected wave BC propagates 
with a velocity equal to 86 percent of the linear-elastic case. 
However, it is shown in Section 5.6.2 that apparently equivalent 
test conditions produced large variations in wave velocity. 

Following the passage of wave BC, the velocities of the gage 
marks in general reduce to values quite small relative to the 
impact velocity and at later times a trend toward the impact 
velocity is visible. As seen in the diagram resulting from 
the graphical wave plotting technique for the linear-elastic 
case (Figure 6l)* the wave action becomes quite complicated 
above the wave BC. Due to the limitations imposed on resolution 
of the motion by the camera framing rates, this complex wave 
motion cannot be traced above the line BC. 

Transverse Waves 

In addition to the longitudinal motion discussed above, a 
lateral motion of the specimen is introduced by the action of 
the impact wheel hook on the specimen. 

The x-dash line in Figure 69 is the locus of the leading edge 
of the resulting lateral wave that propagates along the speci- 
men with an initial velocity that is considerably less than 
the previously determined velocity of propagation of longitudinal 
waves. Extrapolation of this locus back to the left boundary 
gives t =0.9 milliseconds, the same value obtained by extra- 
polation of wave line AB. 

It can be seen in Figure 69 that the locus of the transverse 
wave front from A to F (in the bridle portion of the specimen) 
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is not a straight line, indicating that the transverse wave velo- 
city is not a constant. The reason for the variation in transverse 
wave velocity can be determined by inspection of the graphical 
wave tracing diagram for the linear-elastic case. Figure 6l , where 
it can be seen that there exist reflected waves OH, HI, IK and KF 
that generate variations in the bridle webbing force which in turn 
produce variations in the velocity of propagation of the transverse 
wave. 


The Lagrangian transverse wave velocity 
by: 


(Reference 32) is given 


( 65 ) 


where F is the tensile force in the webbing, p Q is the unstrained 
lineal mass density and e is the strain in the material ahead o_ 
the leading edge of the transverse wave. The velocity relative 
to the laboratory is given by 

=t T = (1 + «) c T ♦ V (- 6 ) 

where V is the longitudinal velocity of the material in which 
the transverse wave is propagating. 


An attempt is made here to explain the curvature in the locus of 
the transverse wave in the bridle by performing a computation 
based on best available estimates of the governing parameters. 

a) The locus of the end of the bridle is estimated 
by drawing a line (GF) parallel to the locus 
of the leftmost gage mark and passing through 
the intersection of the wave AB wi'*h the locus 
(GL) of the bridle end prior to the arrival 
of the wave. 
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b) The actual reflected waves Gil, HI, IK and KF 

are approximated in Figure 69 by using the slope 
of the initial wave AB in their construction. 

c ) The longitudinal strain levels In the material 
into which the transverse wave is propagating 
are computed from the impact velocity, the ob- 
served wave velocity and the longitudinal strain 
wave amplitudes shown in Figure for the linear- 
elastic case. 


S AB 

Vo 

c 

275 

3790 

.0726 

S DE 

.6667 

£ AD 

.0484 

£ EJ 

.3333 

£ AD 

.0242 

S JF 


C AD 

.0323 


d) The velocity of the webbing at any point in the test 
specimen is given by 

v <*> = V o + I[ °< £ R- e L> \ (67) 

where e_ and are the strain levels on the left and 
right side of the ith strain wave respectively and the 
summation is over all strain waves between the left 
(impact) end and the point in question. The webbing 
velocities for each segment of the AF portion of the 
transverse wave locus are therefore. 
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V AD * V o ’ - 27 5 «/»<=« 

V DE * V o + ( ° >AB (0 - 3333) <‘AD J 

- -275 + (-3790) (0.333X0.0726) = -367 ft/sec. 

V o = -275 ft/sec 

V o + (c )ab (-.1111) U AL ) 

-275 + (-3790) (-.1111) (.0726) = -245 ft/sec 

e) The effective stiffness of the bridle is taken to be 
that value of bridle stiffness which gives the observed 
AB longitudinal wave velocity in a linear-elastic 
webbing 



Solving for k, 

k = pc 2 ( 69 ) 

For the bridle, p = .00115 slugs/ft so that 
k = (.00115) (3790 ) 2 = 16,500 lb 

f) With the addition of the linear expression for force 

F = ke ( 70 ) 

all necessary information for the computation of c , the 

transverse wave velocity with respect to the laboratory, 
is available. Substituting into Equations ( 6 s ) and 
(6 6 ) there is obtained 
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(<U > 
t ad 

If 

00 

0 

ft/sec 

<v 

= 485 ft/sec 

T DE 



(c* > 

^T EJ 

= 320 

ft/sec 

( C ^JF 

= 446 

ft/sec 


g) Using these values of c , the locus of the transverse 

T 

wave can be constructed in stepwise fashion. The result 
is shown in Figure 70 as the dash line ADEJF. The 
actual locus is shown as x's. 

The close proximity of the two curves implies that for the type 
of webbing used in the bridle and for the stress levels existent 
in these tests. Equation (£g ') gives a good prediction of the 
transverse wave velocity and the noted curvature of the transverse 
wave locus in the bridle is to be expected. 

At Point F (Figure 69 )* there occurs an increase in the 

transverse wave propagation velocity. It is apparent from the 
diagram that three contributory events occur almost simultaneously 
at Point F. 

a) The transverse wave passes from the bridle into 

the specimen trunk. The wave, therefore, moves into 
material with one-half of the lineal mass density of the 
bridle. 
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Time, milliseconds 


MRIHMP 



Distance From Deadman, in. 

Fig. 70* Theoretical and Experimental Loci of the 

Transverse Wave in the Bridle, Test Type D5. 
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b) The positive strain wavt 3C arrives. 

c) Ti.e positive c trail. XF -.rriv-r-c. 

As can be seen iron Equations -cj and 7C, all three of these 
events have the effect of Increasing the transverse wave propagation 
velocity. The multiplicity of wave fronts makes the application 
of the methods of calculation of the preceding section impractical 
In the trunk of the specimen. 


Further Observations 

The two breaks are indicated by circles in Figure 69* The 
first failure occurs at the right end of the trunk section 
before the arrival of the transverse wave. The location of 
the second break with respect to the first implies that the 
second failure cannot be affected by the first without the 
imposition of unrealistic wave propagation velocities. The 
possibility of such multiple independent failures was indicated 
by the results of the approximate linear-elastic analysis in 
Section 5.^. 

p.6.2 Variations In Longitudinal Wave Velocity and Strain 

The films of similar tests on six additional specimens with 
the same trunk material as Test D5-1 were analyzed to determine 
the apparent propagation velocity of the first wave through the 
trunk and the average strain in the trunk after the passage of 
the first wave. The epace-time diagrams obtained from the 
reduction of the respective film data are shown in Figures 71 
to 76. The dashed lines in these diagrams are the best fit 
straight ^ines indicating the apparent wave motion and the 
motion of the gage marks. The apparent wave velocity, the 
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Fig. 71 . Space-Time Diagram, Test D 3 - 
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Fig. 74. Space-Time Diagram, Test D4- 
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average strain after the wave, and the average particle velocity 
after the wave are displayed in Table n. The entries are 
arranged in order of increasing average strain. 

Table 11 shows large variations in apparent wave velocity and 
strain that do pot appear to have any simple correlation with 
the only supposed variations between the tests, the impact 
velocity and the matching of the webbings in the trunk, Rwever, 
there does appear to he a significant level of correlation, cetween 
apparent wave velocity and strain; apparent wave velocity pends to 
increase as the strain decreases. 


This trend is to be expected because according to References 36 

and 37 the particle velocity v^ attained at station i in the 


specimen in given by the expression 


v 


i 


\ c(e)de 

,J 0 


( 71 ) 


where c is the longitudinal strain at station I and where c(e) 

Is the wave propagation velocity at the strain level e. 

The velocity c is given by 

( \ f 1 

c < £ > = (p 0 3a 1 ( 72) 

where F is the tensile force in the specimen and p Q is the lineal 
mass density in the unstrained condition. It can be seen from 
Equation 71 that for approximately equal values of particle 
velocity v^ large Increases in wave propagation velocity must be 
accompanied by an appropriate decrease In the strain level attained 
behind the wave. 
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Table n . Comparison of Wave Velocities, 
Particle Velocities and Strains 


Test No. 

Impact 
Velocity 
(ft/sec ) 

Wave 

Velocity 
(ft/sec ) 

Particle 
Veloc ity 
(ft/sec ) 

• “i 

Strain 

' 

Particle 
Velocity, 
Eqn. 71 , ! 
Static 
(ft/sec ) 

D3-3 

-300 

4700 

-310 

.069 

-250 

D4-6 

i 

-285 

5550 

-317 

.070 

-252 | 

D3-1 

-284 

| 3040 

-234 

.080 

j 

i ’ 293 j 

D6-1 

-277 

3850 

-272 

.090 

-340 1 

D5-1 

! 

- 275 i 

| 

3790 

-294 

.091 

-342 | 

| 1 

D4-5 

-277 

3230 

| 

-313 

) 

1 . 093 

[ 

j ~ 355 ; 

D3-2 

-284 

1 2550 

| -268 

.123 

-550 i 
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References 22 and 25 show that the dynamic load-strain curves, 
for some typical nylon yarns for the case of monotonically 
increasing strain at constant strain rate, approach a linear 
type of relationship as the rate of loading Increases. Also, 
the slopes of the dynamic load-strain curves tend to increase 
with rate of loading. If the dynamic load-strain curves fol- 
lowed by the materials in these tests during the passage of the 
first wave are linear, then equation i71 ) becomes 


Vi = cei 


(73) 


or 


c = 


Vi 


(74 ) 


i.e., the product of the wave velocity and the final strain 
must be equal to the particle velocity in all tests. 


Figure 77 is a plot of the left side of Equation 74 as a 
function of the quotient on the right side; one point is obtained 
for each test. If in fact the specimen trunk material did per- 
form in a 3 inear fashion during passage of the first wave then 
all of the points should fall on a line of unit slope passing 
through the origin. The scatter that is obtained .is of the 
order of magnitude of the inherent errors in the determination 
of velocities by graphical differentiation of the data. How- 
ever, there does appear to be a trend above the unit slope which 
could possibly be caused by a consistent small nonlinearity in 
the performance of the material. 


It is now demonstrated that the performance of the trunk 
material is not governed by the static load-strain curve. This 
is done by using Equation 71 to compute the particle velocity 
using the static load-strain curve to compute c(e) and the 
value of ^ attained behind the first wave. Since the static 
load-strain curve is concave upwards, shock wave effects must 
be taken into account as shown in References 35 and 37. 
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0 2000 4000 6000 

V 0 / € , ft/sec 


Fig. 77 . Observed Longitudinal Wave Velocity 
vs. Wave Velocity Calculated on the 
Basis of a Linear Load-Strain Curve 
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Reference 37 gives an excellent set of examples of such calcu- 
lations. The resulting values of particle velocity are given 
in the last column of Table Hand there appears to be no cor- 
relation with the actual particle velocity as recorded in the 
fourth column of the same table. The implication is that the 
static load-strain curve is not in itself a valid predictor 
of the dynamic performance of nylon webbing. 

As yet unresolved is the question of the source of the mutually 
consistent variations in wave velocity and strain. At the time 
of the writing of this report, it is believed that, even though 
the impact velocities were approximately equal, the rate at 
which the end of the specimen was accelerated to this velocity 
was different in each test due to variations in the construction 
of the impact loop, the initial shape of the loop and the manner 
in which it was engaged by the hook. If the dynamic load-strain 
curve of the specimen trunk material is a sufficiently strong 
function of the rate of loading, variations in apparent wave 
velocity and strain such as found in these test". are possible. 
Attempts to correlate wave velocity and strain with test en- 
vironmental conditions and the amplitude of the transverse wave 
were fruitless. 

5.6.3 Estimation Of The Dynamic Load-Strain Relationship 

The velocities of the longitudinal and transverse waves and 
the longitudinal strain can be used to calculate the force 
in a specimen trunk in two regions of the space-time diagrams. 
When divided by the number of webbings in the trunk and plotted 
against the longitudinal strain in the respective region, the 
result is three points (one is the origin) on the dynamic load- 
strain curve followed by the material of each trunk in each 
test. Several assumptions are involved in the calculations. 
These assumptions are presented at the appropriate place in 
the following description of the procedure. 
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If It Is assumed that the longitudinal wave shape remains con- 
stant as the first wave progresses across the specimen trunk, 
then, using the principle that the net force acting or the 
material particles containing the wave is equal to the rate 
of change of momentum introduced by the waves gives the fol- 
lowing expression for the force in the specimen behind the 
first wave at station i. 

=-Po c iVj[ (fo) 

Table 12 shows the quantities used in the calculation of force 
by this method and the result. The second and sixth columns 

give a single point on each load-strain curve. 

Equations 65 and 66 can be solved for the force in the webbing 
to obtain 

P = v c tT - V, 2 <T6) 

This equation is used to compute the force in the webbing at 
the regions in space-time where the transverse wave locus inter- 
sects relatively smooth gage mark loci. The ingredients in the 
calculations and the results are displayed in Table 13 . 

The last two columns give an additional single point on each 
load-strain curve. 

The resulting partial dynamic load-strain curves are shown in 
Figure 78 Two significant observations are that the dynamic 
load-strain curves in some cases lie below the static load- 
strain curve (shown as a dashed line) and that the material 
does not always follow the same curve in these tests. It is 
felt that the inherent scatter in the data causes the apparent 
variations in the sign of the curvature of the load-strain 
curves and does not permit the detailed inspection required to 
make conclusions regarding the degree of linearity of the 
dynamic load-strain curves. 
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Table 12 Determination of Points on the Dynamic Load -"train Curves 

Using the First Longitudinal Wave Across the Specimen Trunks 


Test No. 

Average 

Strain 

Lineal 

Density 

(sl/ft) 

Wave Velocity 
(ft/sec ) 

1 

Average 
Par f ic le 
Velocity 
( ft/sec ) 

Force/ 

Webbing 

d v ) 

D3-1 

.080 

.00118 

3040 

-234 

409 

D3-2 

.123 

.00113 

2330 

- 2 f ’8 

393 

D3-3 

.069 

.003 13 

4700 

-310 

8 38 

D4-5 

.093 

.00115 

3230 

-313 

381 

D4- r 

.070 

.00115 

3550 

- -07 

1010 

D5-1 

.091 

. 000374 

3790 

-234 

640 

DG-1 

.090 

.00115 

7850 

-272 

^02 


Table 13 Determination of Points on the Dynamic Load-Strain 
Curves Using Transverse 'Wave Velocity 


Test No. 

Gage 

Marks 

Used 

Wave 

Velocity 
With 
Respect 
to Labo- 
ratory 
( ft/ sec ) 

Particle 
Velocity 
(it/ sec ) 

Strain 

Force / 
Webbing 
(lb) 

D3-1 

D3-2 

* 

1,2 

1486 

-109 

0-197 

1220 

D3-3 

2,3 

1625 

-106 

.186 

1452 

D4-5 

2,3 

1690 

-107 

0.145 

1621 

D4-6 

2,3 

2170 

-117 

O.loO 

2592 

D5-1 

2,3 

1880 

- 99 

0.184 

I 885 

D5-1 

3,4 

1670 

-155 

0.218 

1572 


* Transverse wave locus does not Intersect gage mark loci 
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78 Dynamic Load-Strain Curves for the Specimen Trunk 
Material as Determined From Wave Velocities 
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It should be emphasized that the load -strain curves displayed 
In Figure 78 are not for the usual case of monotonically in- 
creasing strain at constant strain rate. In fact, between the 
states Indicated by the two points on each curve, the space-time 
diagrams indicate that for a significant interval the strain 
remains essentially constant. 

5.7 SUMMARY OF PILOT CHUTE RISER DYNAMICS STUDY 

It is shown that many aspects of the pilot chute riser test da 1 :' 1 
can be explained with a linear-elastic approximate analysis. 

It is clear from this part of the study that wave phenomena 
must be considered in the structural analysis of parachute 
risers. That is, the mathematical model of the riser must take 
into account that the riser is a continuum, rather than a small 
number of lumped springs and masses. It is also shown that 
multiple failures and failures in free lengths are distinct 
possibilities in destructive tests and that the equation of 
state used in a riser dynamic analysis must take hysteresis 
effects into account. Comparison with experimental data shows 
that the linear-elastic approximation Is not adequate for riser 
structural analysis. 

A mathematical model which includes the inertia effects of the 
Impact wheel, motion of the webbing around the Impact wheel and 
the nonlinear static load-strain curve as the equation of state 
produces improved results. This mathematical model is analyzed 
using an existing finite-difference method computer program 
(AFTON) which is extensively modified to be compatible with the 
boundary conditions of the tests and to accept the nonlinear 
static load-strain curve of nylon webbing as an equation of 
state. Further refinement of the mathematical model is found 
to be necessary. 
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A space-time diagram of one riser test is analyzed in detail 
with respect to the motion of the longitudinal and transverse 
waves and several other diagrams are analyzed in less depth. 

It is shown that the dynamic load-strain curves followed by 
the nylon webbing are approximately linear but small variations 
from linearity produce effects which cannot be neglected. The 
longitudinal wave speed appears to be sensitive to some as yet 
undiscovered variant in the tests, perhaps the nature of the 
impact loop. 
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SECTION 6.0 


MEASUREMENT OF PRESSURES, LOADS, AND STRESSES 

The types of measurements that are required to support the 
further development of the Internal load prediction methods 
being developed for Apollo type parachutes are discussed In 
this section. These measurements are needed to provide Input 
data for the Internal loads analysis and to confirm the accuracy 
of the analysis. 

During the current Apollo analysis study, the CANO stress 
analysis program has been refined Into a powerful analytical 
tool. Section 6.1 discusses the CANO program Input data 
requirements and shows how the experimental measurements to be 
made under the proposed test program will both provide these 
data and serve as a check on the validity of the method. 

6.1 STRUCTURAL ANALYSIS DATA REQUIREMENTS 

The analysis for determining the Internal load distribution In 
a parachute under known external loading Is given In Section 3. 

Input data to CANO consists of parachute construction geometry, 
material load-strain curves, canopy differential pressure 
distribution, and applied riser load. Output includes loading 
and strains in radial tapes, horizontal members (sails or ribbons), 
suspension lines, vent lines, vent band, skirt band and reefing 
lines. In addition, a complete geometric definition of the 
parachute shape is obtained. 

The CANO computer program Is limited to symmetrical shapes, and it 
is assumed that the applied riser load Is reacted by differential 
pressure acting normal to the local canopy surface. 
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For a given oarachute construction there is a unique relationship 
between the internal stresses* the aoolied loading, and the 
geometry of the deformed canopy elements. This relationship is 
represented in the diagram of Figure 79a. Program CANO is arranged 
to proceed in the direction shown by the arrows in the diagram. 

The external loading, consisting of the riser load and a nondimen- 
slonal pressure distribution curve, and the parachute construction 
details comprise the input data. The internal loading, consisting 
of a load for each element of the structural model, is the primary 
output. Definition of the canopy equilibrium shape and strains 
is a by-product of these computations. 

The information in one block of Figure 79a is sufficient for deter- 
mination of the other two. However, the accuracy of measurement 
that would be required to infer internal and external loading 
from the canopy geometry is beyond the practical limits of 
photographic analysis. 

In the test program described in this report, it is planned that 
all of the parameters shown in the three blocks of Figure 79a 
will be measured. The redundancy thereby obtained will provide 
a check on the validity of the assumptions made in derivation 
of the analysis. 

The external loading to be measured consists of the riser load 
and the differential pressure acting across the canopy surface 
The riser force can be measured accurately by load links, tut 
measurements of the differential pressure presents many difficul- 
ties. The practical limit for the number of pressure transducers 
is a test is about ten. In a ringsail parachute the pressure 
varies locally between the leading and trailing edge of each 
sail in addition to the general variation from skirt to vent. 

The redundant measurements of stress, strain, and shape will be 
valuable in deriving, from the limited pressure measurements, a 
pressure distribution curve which will give an accurate stress 
distribution for use in structural analysis. 
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Fig. 79. CANO Program Input/Output Relationships 
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If the external load block of Figure79a is considered as two 
separate data groups, pressure distribution and riser load, the 
diagram shown in Figure79b is obtained. 

In Figure 79^ any two blocks provide sufficient information for 
computation of the other two. The riser load measurement is the 
simplest and most reliable of the four. If it Is assumed that 
this value is always measured accurately, the requirements for 
the other three measurements are less suringent For the pres- 
sure differential and stress, only the shapes of the distribution 
cures are required. For the shape data, only the relative shape 
of the radial tape profile is necessary. This Is helpful because 
absolute values for length measurements from photographs are 
the most difficult to obtain. 

Since the magnitude as well as the distribution of the differential 
pressure will be measured, additional redundancy exists. The 
CANO program requires only the shape of the pressure distribution 
curve, the magnitude being automatically adjusted to balance the 
riser load. Comparison of the pressures computed in this manner 
with the measured values will give a check on the analysis, in 
particular the assumption that skin friction and drag of fluttering 
sails is negligible. 

The simultaneous measurement of loads and strains will provide 
Information on the dynamic and repeated loading behavior of 
the parachute materials. Although the CANO program is not limited 
to static load-strain curves, use of dynamic data is not considered 
worthwhile at present in view of the uncertainties in other input 
data. 
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6.2 SUMMARY OF PRIOR WORK ON CANOPY PRESSURE MEASUREMENT 

The first measurements of pressure distribution in a parachute 

2 

canopy were made by Jones in 1923 for the Advisory Committee 
for Aeronautics of Great Britain. A 4-foot diameter wind 
tunnel at the National Physical Laboratory was used for the 
measurements, with sting-mounted models constructed of wood, 
copper sheet, and silk stretched over a metal ribbed framework. 

Models were approximately 8 in. in diameter and simulated a 
canopy shape developed by Taylor for the Advisory Committee 
for Aeronautics. Steady state pressures were measured by insert- 
ing two lengths of hypodermic tubing in each model surface, one 
flush with the interior surface and the other with the exterior 
surface One end of the tube was plugged and holes were drilled 
at each pcint where pressure was to be measured. The measurements 
were recorded on a manometer, by uncovering one hole at a time, 
at velocities from 30 to 50 ft/sec and angles of attack from 0 
to 10 deg. 

Modern efforts at measuring pressure distribution in a steady 

state canopy have been carried out by Heinrich and by Babish 

and Hunter^. Heinrich’s work also is described in the Air Force 

15 

Parachute Handbook . Stainless steel, 2.5-inch diameter models 

of flat circular ribbon and guide surface canopies were oi.ing 

mounted and wind tunnel tested at Mach numbers between 0.6 and 

5 5 

1.25 over a Reynolds number range from 7 x 10 v to 9.7 x 10 . 
Measurements of internal and external pressure distribution were 
made, and shadowgraph pictures were taken in an effort to aid 
visualization of the flow field through ribbon canopies. External 
pressure coefficients in the subsonic flow regime (below Mach 
number 0.8) were found to be relatively insensitibe to Mach 
number. Internal pressure was found to vary little from free 
stream stagnation pressure at all points in the canopy. 
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The work carried out by Babish and Hunter was a cleverly conceived 
effort to arrive at three-dimensional canopy pressure distribution 
by pressure and flow field measurements on two-dimenslonsal rigid 
models. These models represented the profile shapes of the 
inflated drag-producing surfaces of various types of parachutes 
as determined from photographic records of descending canopies. 

The various parachute types represented included flat circular, 
flat circular ribbon, ringslot and ringsail, and others. The 
models were tested in the Air Force Institute of Technology Smoke 
Tunnel and were constructed so as to completely span the test 
section. Test section velocity was 29 ft/sec # which provided 
ideal smoke streamline flow. By measuring tht steam tube thick- 
ness and utilizing a theoretically developed relationship between 
streamtube thickness and pressure, the pressure coefficients 
could be obtained in regions of the flow where the streamlines 
were well defined. To obtain measurements in regions where the 
streamlines were not well defined, i.e., the stagnation and 
wake regions, an instrumented semicircular concave cylinder was 
tested. By combining the above two types of measurements, the 
pressure distribution over the entire surface of the two-dimensional 
models could be obtained. These pressure distribution data, 
reduced to appropriate force coefficients, were compared with 
previously established force coefficient data for three-dimensional 
models. By use of this method and appropriate assumptions, it was 
possible to adjust the two-dimensional data to yield the pressure 
distribution about an actual canopy. Difficulties were encountered, 
however, in that the streamlines close to the body entered the 
boundary layer flow, causing errors of + 20 percent in streamtube 
thickness, and resulting in a corresponding lack of confidence in 
the pressure distribution. The primary value of this study is 
considered to be the excellent smoke streamline photograpns and 
the insight which they afford into the steady state flow about the 
canopy. 


188 


NVR 6432 



The reports reviewed above have been concerned with the pressure 
distribution and flow field about a fully Inflated, steady state 
canopy. Although these studies are Important because of the 
Insight they give into the operation of a steadily descending 
parachute, the object of primary interest in the current study 
is the pressure distribution in an Inflating canopy, since it is 
this phase of the process which Imposes maximum stresses on the 
canopy . 

15 

Hefnrich obtained measurements of canopy pressure distribution 

on a series of fixed models representing the shape of the canopy 

at various stages of Inflation. The results of this study were 

4o 

used in a subsequent paper by Heinrich and Jamison on predicting 
canopy stresses during inflation. The tests were conducted on 
seven synthesized shapes at a low subsonic speed, and the pres- 
sure distributions associated with each shape were presented in 
the report. Although these tests may indicate the proper trends 
for demonstration of stress analysis techniques, it is doubtful 

that they satisfy '-he need for pressure measurements under dynamic 
conditions . 

The only Investigations made to date wherein pressure distribution 

20 21 4l 42 

was measured in an inflating canopy are those of Melzig * » * 

ft 

of the Deutsche Forschungsanstalt fur Luft and Raumfahrt (DFLR), 
Braunschweig, Germany. Of the four published works by Melzig, 
three involve actual drop tests and one concerns a wind tunnel 
test program under infinite mass conditions. The latter will be 
discussed first. 

The wind tunnel investigation was made in a 9 x 12 foot low speed 
open- jet wind tunnel at velocities between 70 and 160 ft/sec. The 
test models were of the following designs: 

(1) Solid cloth, circular flat - 28 gores, D Q « 53.5 in. 

(2) Solid cloth, 10# flat extended - 30 gores, D Q * 62 in. 

(3) Ringslot, geometric porosity 16# - 24 gores, « 53.5 in. 
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The models were mounted In the tunnel in a stretched-out position 
with leather clamps for retention. After the tunnel was brought 
up to speed, the models were suddenly released. Pressure distri- 
bution (internal, external and differential) was measured in the 
canopy by pressure transducers attached at four specific locations. 
The pressure transducers were developed at DPLR, and were 1.2 in. 
in diameter, 0.35 in. thick, and weighed 0.2 oz each. The sensors 
had a capacity of + 0.5 psi, and were compensated for temperature 
and for acceleration to 200 g's. The filling process was photo- 
graphed from one side of the tunnel at a camera speed of 100 frames 
per second and individual frames were synchronized with the 
pressure recordings. Four identical tests of each canopy type 
were made at each selected test velocity. It was found that the 
reproducibility of the four measurements was reasonably good. 

A canopy shape analysis was made based on the photographs obtained. 
Various phases of filling were .defined by assigning an "ideal 
photographic shape" made up of appropriate bodies of revolution 
as a function of time were presented in graphical form for all 
four canopy types. 

The remaining references on Melzlg's work concern canopy pressure 
distributions and inflation profiles obtained during free flight 
drop tests. Reference 41 reports the measurements made during a 
series of nine drop tests of a T-10 parachute at speeds from 80 
to 140 knots with a 200 lb weight attached. Differential pressure 
transducers identical to the type described above for the wind 
tunnel program were mounted at four locations in the canopy between 
the skirt and the vent. The pressure showed considerably more 
oscillation and lack of a characteristics pattern than had been 
the case in the wind tunnel testing presented in Reference 11. 
Reference 42 describes results of 18 drop tests on 28-foot dia- 
meter circular flat and hemispherical parachutes at an initial 
drop velocity of 118 knots. Four differential pressure transducers 
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were again placed in the test canopies and their output recorded 
on an oscillograph in the drop dummy. Ground photographs of 
the inflating canopies were made at a speed of 100 frames per 
second by a 35 nun camera, and also by a camera mounted aboard 
the dummy. The pressure traces again showed severe oscillations 
in the initial (flutter) phase of the opening process. A 
second phase was entered in which the negative pressure peaks 
disappeared and the oscillation frequency decreased. Following 
this was a third phase in which a positive steady pressure 
appeared at all points in the canopy. Melzlg correlated the 
pressure trace behavior with the canopy opening history and 
load traces. 

21 

The most recent work by Melzig involved drop tests of four 
parachute types j circular flat (C-9), 10 percent extended 
skirt (T-10), ringslot (two designs) and circular flat ribbon 
(two designs). The first ringslot had a nominal diameter of 
32 feet, 32 gores, 10 cloth Btrips, and 15 percent geometric 
porosity. The second had the same dimensions but 11 instead 
of 10 cloth strips. The circular flat ribbon parachutes had 
a nominal diameter of 32 feet, 34 gores, and a geometric porosity 
of 15 percent. One had 83 and the other 76 concentric ribbons. 
Both onboard and ground cameras were used to photograph the 
opening process at 100 frames per second. Four DFLR differential 
pressure transducers were equally spaced in the canopy from skirt 
4 to vent. Riser force, vehicle dynamic pressure, and the four 
differential pressures were telemetered to ground receiving 
equipment. The vehicles were dropped from an altitude of 
1000 feet at a velocity of 110 knots. 
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In general, the pressure histories show the same overall trends 
for all four types of parachutes. There are rapid fluctuations 
between positive and negative peaks during the flutter phase. 

The fluttering ends relatively soon at the vent, at which point 
canopy filling begins, traveling progressively toward the skirt 
as positive pressure is gradually established over the whole 
canopy. 

In conclusion, Melzlg states that because of the dynamics of the 
inflation process only qualitative information on canopy stresses 
can be obtained from the detailed pressure distribution measure- 
ments, and that idealization cf the shape during the flutter 
phase can yield false results. Additionally he states that be- 
fore the pressure measurements during the flutter period can be 
usefully applied, stress analysis techniques must be developed 
which account for cloth dynamics. 

Despite Melzlg 's emphasis on cloth dynamics during the flutter 
phase, experience with Apollo parachutes has shown that canopy 
stress analysis methods which consider only filled sails have 
accurately predicted test failures. Therefore, special efforts 
to measure canopy pressure, strain, load and shape during the 
flutter phase are not warranted. 

6.3 PRESSURE MEASUREMENT 

Several types of pressure transducers have been proposed by 
various manufacturers for parachute application. Some manu- 
facturers have also proposed complete pressure measurement 
systems including telemetry and recording instruments. However, 
the tests by Melzlg, described in Section 6 . 2 » are the only 
ones in which inflight pressure measurements were made. The 
Melzlg type transducer is the only one which can be considered 
flight qualified. The following paragraphs describe the various 
pressure transducers and measurement systems which have been 
proposed . 
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SchJeldahl System 

Schjeldahl Company, Northfiexd, Minnesota , has proposed a dif- 
ferential pressure monitoring system (Pigure80a) which does not 
require wires in the parachute’s risers. The basic sensor is 
similar in operation to that developed by Melzig at DFLR 
except it offers improved environmental performance. However, 
packaged within the same unit are the sensor, electronics for 
& miniature PM transmitter, four small mercury batteries and 
signal conditioning electronics. Each unit is 2-1/2 in. in 
diameter, 1 in. in thickness, weighs approximately 1.75 oz and 
requires no external wires for operation. The miniature trans- 
mitter sends the sensor data to either a receiver-recorder or a 
repeater (which transmits the Information to the ground) located 
in the payload. Lanyard switches are proposed for activating 
the units upon drop of the test vehicle. The system is limited 
to a total of seven transducers by telemetry charnel availability. 

RdP West System - A system using a different type of sensing 

lijl 

element has been proposed by RdP West, Westminster, California. 

Their system (Figure 80b) uses a variable, stretched diaphragm 
capacitance transducer and solid state signal conditioning elec- 
tronics to provide a DC output pressure transducer system. The 
sensors are 1.12 in. in diameter, 1 in. in thickness, weigh 1 oz 
and exhibit accuracies comparable with other sensing devices. 

The master signal conditioner, located at the skirt of the para- 
chute, is 1.25 x 1.25 x 2.00 Inches am! weighs 4 oz. The con- 
ditioner output is fed into small voltage controlled oscillators 
and a mixer amplifier. This unit, also located at the skirt, has 
an approximate volume of 4.5 in^ and weighs 4 oz. The wire carrying 
the mixed PM signal is run down a suspension line and connected to 
the tape recorder and/or transmitter. Three power supply lines 
are the only additional wires needed on the suspension lines. 
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Fig. 80 • Schjeldahl and RdF West Pressure Measurement Systems 


194 


NVR-6432 




Manning Instruments Transducer - Manning Instruments Division 
has proposed a transducer of novel construction for measuring 
differential pressure in a high acceleration environment. This 
transducer Is shown schematically in Figure 8l. It consists of 
twin diaphragms with the external surfaces exposed to the higher 
pressure level. The inner surfaces are exposed to the lower 
pressure level. T*o silicon strain gages are bonded to the 
external surfaces of each diaphragm. The gages on each surface 
are placed (Figure 8l)so that gages A and C are strained in 
compression and gages B and D in tension when the transducer 
is loaded normally. The transducer is self-compensated for 
acceleration in that acceleration causes both diaphragms to 
deflect In the same direction, thereby balancing the bridge by 
straining gages A and B in the same direction and gages C and 
D in the same direction. This transducer is available in dif- 
ferential pressure ranges between 5 and 100 psi and weighs 
about 0.3 ox. 

Sensotec Transducer - The Sensotec absolute pressure transducer 
is a miniaturized capsule (.25 in. diameter x .03 in. thick), 
one side of which is the active diaphragm. The diaphragm has 
a pressure sensitive area of .023 in 2 and utilizes semiconductor 
8 train gages bonded to the inner surface. The differential 
pressure model has the same diameter, but is about 0.2 in. thick. 
The smallest differential pressure range currently available 
is 0-2 psi. Although no acceleration compensation is built into 
the transducer, Sensotec claims errors due to acceleration at 
200 g would be less than 1 percent. A temperature compensation 
network can be built into the mounting pad. 
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Solid State Strain Cages 
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(a) Sensor Schematic 



(b) Wiring Diagram 


Fig. 81 • Manning Instruments Pressure Transducer 
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6.4 STRAIN TRANSDUCERS 

Strain Transducer Based on Inductive Circuit - During the present 
study, a concept was formulated for a canopy strain measurement 
device based on the Inductive effect between two coils located 
near each other (Figure 82). The two colls are orthogonal, and 
an output null point can be realized. The coils can be mounted 
on the canopy cloth so that they move relative to each other as 
the cloth stretches. The output signal Is dependent on the 
physical displacement of the coils, which Is In turn a measure 
of the strain In the fabric between the attachment points of the 
two colls. 

Hie electrical characteristic upon which this transducer Is 
based Is commonly used in resolvers and similar electronic 
equipment. However, a laboratory test program is required to 
develop a device of acceptable accuracy, size and weight. 

NASA Langley has performed strain measurements by bonding 
so-called post yield strain gages to the parachute tape members. 
Strain gage manufacturers claim that these gages which are made 
of fully annealed Constantan, can be used to measure 10-15 per- 
cent. Tests at NASA Langley have indicated good repeatability 
in measuring 10 percent strains during short duration dynamic 
load pulses. The Langley tests have shown 1 to 2 percent residual 
In the fabric and not due to a gage shift. 
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Pig. 82 . Variable Inductance Strain Transducer 




Non-electrleal Strain Sensor 


A non-electrical maximum strain indicator for high elongation 

45 

materials has been developed by R. L. Hanes at Northrop Ventura. 
Figure 83 shows the principle of operation of this device. 

Laboratory tests conducted with these strain Indicators attached 
to nylon parachute fabric indicate the accuracy of the instru- 
ment to be within one-percent of the gage length. This accuracy 
is satisfactory for the measurement of loads or strains In para- 
chute fabrics. A number of these gages have been used in Parawing 
flight tests* but results have not been published. Laboratory 
tests under dynamic loading conditions have not been made. 

6.5 LOAD TRANSDUCERS 

Consideration was also given to a strain gage type transducer 
capable of measuring the load in the radial s of rlngsall type 
parachutes. This device is sketched in Figure 84. The transducer 
is intended to avoid stress concentrations in the material to 
which it is attached. This ts accomplished by attaching identical 
devices to both sides of the ribbon or radial* and by incorporating 
stiff rubber mounts for holding the metal plates to which the 
strain gages are attached. The rubber stiffness is selected to 
be somewhat greater than that of the ribbon or radial to which 
it is bonded. This allows the entire load to be transferred to the 
strain plates without causing stress concentrations in the ribbon 
material . 

A development program would be required to produce a transducer 
of this type that would be suitable for measuring parachute loads. 

Considerable work has been carried out at NASA Langley on 
development of a miniature load transducer that is applicable to 
tape members of parachute canopies. This device can be sewn to 
the tape member at whatever location is desired. It consists of 
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Elastic tape backing attached 
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Load relieved condition 

Fig . 83 . Nonelectrical Maximum Strain Sensor 
(Reference ^3 ) 


800 


NVR-4432 





201 


NVR-6432 


Pig. 84 • Radial Tape Load Transducer 



a rectangular stainless 3teel plate instrumented with four small 
strain gages of standard design and having end tapes for 
attachment to the radial. Loads up to 180-lb have been measured 
at Langley in wind tunnel and flight tests. 

6.6 PHOTOGRAPHIC INSTRUMENTATION 

Inflight Testing - Accurate profile and plan view shapes of 
inflating parachutes are among the most valuable data to be 
obtained; however, these data also seem to be among the most 
difficult to obtain. Cameras mounted on the test vehicle 
obtain reasonably good plan view pictures of parachute canopies. 
These cameras can be calibrated by photographing a grid of 
known dimensions. The distance from the camera to the canopy 
can be fairly accurately estimated, so that determining true 
dimensions from such photographs is possible to perhaps 5 per- 
cent accuracy. 

Good profile views seem to occur mainly by chance during flight 
drops, primarily because it is difficult for photo chase aircraft 
to maneuver so as to maintain the proper camera position. Ground 
photographs usually show oblique angle views, making it necessary 
to estimate the camera viewing angle and correct the measured 
profile height. Also, in the case of the drogue chute, the 
resolution seldom is high enough to yield a profile shape suitable 
fbr stress analysis purposes. 

It Is beleived that the best way to overcome these problems in 
the proposed test series Is to test at the lowest possible 
altitude. Ground cameras such as the El Centro "Big Eye" have 
focal lengths such that good resolution can be obtained from 
long distance, so that the drop can be made far from the camera. 
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By making the drops at low altitude, the obliqueness of the view 
i3 reduced to a minimum. After a first order correction, the 
resulting profile shapes should be of sufficient accuracy to be 
used for comparing with shapes theoretically predicted by the CANO 
stress analysis program. 

In summary, it is believed that presently available photographic 
recording means (on board cameras and long focal length ground 
cameras) are adequate provided the test drops are performed 
at low altitudes (less than 2500 ft). 

Wind Tunnel Testing - Wind tunnel tests provide the best 
opportunities for obtaining good plan and profile views of the 
inflating parachute. The plan views can be photographed 
accurately using cameras mounted on the support pylon. Profile 
views, are the most difficult to obtain. One problem is that 
the rays of light from the image to be photographed are not 
parallel . If the distance from the parachute to the camera were 
infinite, or if an optics system capable of producing parallel 
rays were available (such as the Schlieren system commonly used in 
wind tunnels), true projections could be obtained. In the 
absence of such equipment, it is necessary to make corrections 
for the parallax error using the known distance from the camera 
to the parachute. 
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6.7 PROGRAM PLAN FOR MEASURING PARACHUTE PRESSURE DISTRIBUTION, 
LOAD, STRAIN AND SHAPE 

This section describes a program for measurement of the parachute 
canopy parameters required for use in the internal loads analysis 
methods described in Section 3.0. It encompasses a three-part 
plan designed to follow a logical sequence of (1) laboratory 
development and testing of suitable instrumentation, (2) wind 
tunnel testing to obtain the required measurements and confirm 
the adequacy of the instrumentation, and (3) limited flight 
tests to obtain further measurements and to correlate wind 
tunnel results. The program is intended as a major step in 
obtaining a complete understanding of canopy stresses. The 
experimental program strongly interacts with and depends on 
the theoretical methods for predicting canopy loads and stresses, 
and one of its most important results would be to confirm the 
accuracy of the theoretical methods. 

6.7.1 Laboratory Development Phase 

A laboratory development phase is considered to be a vital step 
in any program aimed at making significant Improvements in the 
state-of-the-art of parachute testing. Several pressure trans- 
ducers are currently available for making canopy pressure measure- 
ments, and a choice as to which is the most suitable can only be 
made after the proper tests have been carried out. A number of 
new concepts for measurement of strain and load in canopy struc- 
tural members need to be evaluated and improved upon. It is 
* 

not recommended that laboratory evaluation of complete pressure 
measurement systems be carried out at this v?me, since it is 
considered that this effort would be premature prior to careful 
evaluation of the various available transducers. 
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Specifically, it is recommended that the following items he 
included in the laboratory phase to develop and test pressure, 
strain and load measurement devices. 

(1) Design and test a strain transducer based cn 
the inductive principle (Figure 82). 

(2) Design and develop a radial tape load transducer 
based on the principle shown in Figure 84. 

(3) Perform environmental (acceleration, vibration 
and thermal) tests and calibration checks on 
currently available pressure transducers. (This 
work is intended to complement calibration and 
environmental tests presently being performed 

by NASA). 

(4) Perform static and dynamic calibrations of the 
nonelectrical maximum strain indicator illustrated 
in Figure 8 j . 

In the performance of items (l) and (2), the following general 
approach will be followed: 

a) Perform circuit analysis. 

b) Prepare breadboard prototype (s) - using laboratory 
type recording equipment. 

c) Perform laboratory tests. 

d) Make necessary design modifications. 

e) Retest. 

f) Select tentative design. 

g) Perform acceleration and temperature sensitivity tests. 

h) Select final design. 
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i) Select and obtain flight components. 

j) Perform worn temperature checkout of flight 
type system. 

k) Perform acceleration and temperature sensitivity 
tests of flight type system. 

t ) Fabricate a sufficient number of units for wind 
tunnel testing. 

6.7.2 Wind Tunnel Phase 

Preliminary Considerations A wind tunnel program will be 
worthwhile both for proving out the laboratory developed measure- 
ment systems as well as for obtaining useful measurements of 
canopy parameters during the opening process. 

(Consideration was given to the El Centro Whirl Tower as an 
alternative to the wind tunnel for making the canopy measurements 
discussed above. The primary advantages of the Whirl Tower are 
(1) the facility of observation which it affords for tests 
which can be conducted under finite mass conditions and (2) 
the relatively low cost and simplicity of testing. However, 
the Whirl Tower is not appropriate for testing large parachutes 
having reefed stages because the time available for opening is 
not sufficient to allow disreeflng. It would be necessary to 
open the parachute directly to the stage being tested. This 
procedure is unattractive in that it does not represent the 
true opening process, and no further consideration was given 
to the use of the Whirl Tower in the program . ) 

Modeling Considerations Because of the practical limitations 
associated with scaling the parachute opening process, it is 
important in wind tunnel testing to use the largest possible 
scale model. For the Apollo main parachute, the largest wind 
tunnels available should be utilized. The Ames 40 x 80-foot 
Tunnel and the Langley 30 x 60-foot Tunnel both can accommodate 
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large parachute models. The Am<?s Tunnel Is capable of operating 
at dynamic pressures up to 100 psf, while the Langley Tunnel is 
limited to about 60 psf. Because of Its larger size and dynamic 
pressure range, the Ames 40 x 80-foot Tunnel is the most suitable 
for Apollo parachute testing. 

Prior experience has shown that to avoid undesirable tunnel 
blockage effects, the drag area of the parachute model should 
be limited to 15 percent of the test section area. For a test 

p 

section area of 3200 ft and a Cj^ of about 0.8, the maximum 

allowable D 0 is about 28 ft. This means that one-third scale 

models of the main parachute and full scale drogue chutes can 

be tested. Tests of a one-third scale model and a reefed full 

scale model of an early rain , parachute design were successfully 

46 

carried out in th Ames Tunnel in 1963. 

Parachute models canopies should be scaled geometrically to 
duplicate porosity. This means that the effective thickness, 
and width dimensions of the materials should be scaled. 

This presents a problem in the case of the Apollo main parachute 
since the 1.1 ounce cloth is about the lightest nylon parachute 
material for which porosity can be controlled. In order to 
maintain the proper stiffness distribution, therefore, all of the 
materials should retain full scale effective thickness. For 
example, the radial tapes should be of the same weave with scaled 
width, and the strength of the suspension lines should be reduced 
by the first power of the scale factor rather than the square. 

This will give a structure that is too stiff oy the scale factor 
and will make it more difficult to measure strains accurately. 

In the wind tunnel tests, strains can be made comparable to full 
scale by testing at higher dynamic pressures. The tests of 
Reference 46 demonstrated a reasonable similarity in reefed canopy 
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shape between a full scale and a one-third scale model constructed 
of the sane material* as discussed above. It may be assumed 
that the differential pressure distribution of the two models 
were similar since the shapes were similar. (Figure 85) 

Program Outline The objectives of the wind tunnel test program 
are summarized as follows: 

1) To obtain measurements of riser force, pressure, 
strain and load distribution in a full scale 
Apollo drogue chute and a one-third scale main 
parachute model, and to record photographically 
the shape and filling characteristics. 

2) To provide pressure distribution data for use 
in the CANO stress analysis program, and to 
provide strain, load and shape data for com- 
parison with the program's predictions. 

3) To verify the suitability of the selected pressure 
strain and load instrumentation for flight test use. 

The program test schedule is presented in Table 14. 

In the case of the drogue chute, the maximum available tunnel 
dynamic pressure of 100 psf is only half of the maximum flight 
dynamic pressure of 204 psf for high altitude abort. Although 
it would be desirable to simulate the maximum dynamic pressure, 
it should be possible to determine the manner in which dynamic 
pressure influences the measured quantities by testing in the 
range from 20 to 100 psf. Moreover, if the CANO stress analysis 
program predicts results adequately over the available range 
of dynamic pressure, a high degree of confidence is justified 
in its adequacy at design dynamic pressures since the velocity 
is not so high that significant compressibility effects are 
encountered. 
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Reefing Percentage 
Canopy Growth vs. Reefing Percentage 
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Table 14. Wind Tunnel Test Schedule 
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The approach taken for the 28 foot main parachute models is 
to test at three dynamic pressures. (The test schedule is 
based on the assumption that pressure transducers suitable 
for testing at 90 psf dynamic pressure will be available.) 

The highest dynamic pressures listed for the second reefed and 
full open stages are approximately three times the maximum 
flight values for a high altitude abort. 

Testing Techniques and Procedures All models will be deployed 
from deployment bags with the suspension lines stretched out 
horizontally from the pylon mount to the deployment bag position. 
Figure 86 shows the proposed arrangement of the model in the 
Ames 40 x 80 foot tunnel prior to deployment. The deployment 
bag is suspended from the celling and stabilized in the tunnel 
by a pilot chute. 

The celling support lanyard, which is brought jutside the tunnel 
through a ceiling access port, also supports the cutter line 
for the deployment bag tie. Figure 8? (taken from Reference 46) 
is a sketch of the deployment bag and lanyard arrangement. 

Prior to tunnel start, the pilot chute is laid on the tunnel 
floor. As soon as air begins to flow through the tunnel, the 
pilot chute Inflates and assumes its Intended position behind the 
deployment bag. The bag tie cutter is operated manually to 
initiate deployment. The deployment bag and pilot chute are 
extracted from the tunnel through the ceiling access port. Reef- 
ing lines on the model parachutes are cut by remotely initiated 
pyrotechnic reefing line cutters. 
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Fig. 86 . Typical Parachute Installation in 40x80-ft Tunnel 



Instrumentation Instrumentation for all tunnel tests to support 
the structural analysis methods will consist of the following: 

1) Pressure Transducers - Ten transducers will be 
spaced along a radial tape and a mid -gore line. 

2) Strain Transducers - Ten transducers will be 
located in positions corresponding to those of the 
pressure transducers, but on radials and gores 
adjacent to those used for the transducers, so 
that no sail is required to support more than one 
transducer. 

3) Load Transducers - Ten transducers will be located 
on radial and vertical tapes. 

4 ) Riser Force Load Link - An instrumented load 
link of standard design will be used to obtain 
riser force. 

5) Reefing and Vent Line Load Links - Instrumented 
load links will be used to obtain reefing and 
vent line loads. 

6) Photography - The suggested location of motion 
picture cameras Is shown in Figure 86. A 16 mm 
camera operating at 128 frames/sec will be used 
to obtain axial photographs and a second camera 
operating at 128 frames/sec to obtain profile 
photographs of the inflating parachute models. 

A variable speed 70 mm tripod mounted camera will 
be used to obtain high quality profile photographs. 


214 


NVR 6432 



6} (Continued) 

The cameras will be equipped with automatic timers 
to provide film timing mark3, and will be synchronized 
with the oscillograph recordings by a flashxng light 
common to all cameras and wired through the oscillograph. 
Camera calibration will be accomplished by photo- 
graphing a grid placed at the tunnel centerline 
having overall dimensions equal to the parachute 
profile maximum dimensions. Error in measuring 
parachute profile diameter from the photographs 
will be less than 5 percent after correcting for 
parallax. 

Output of the strain and load transducers and the riser, vent 
and reefing line load links will be recorded on oscillographs 
or Sanborn recorders. 

Data Reduction - Scaled illustrations at selected times during 
the opening process will be prepared for each parachute tested. 
All pressure, strain and load data (including riser force and 
vent .line force) will be smoothed and plotted as functions of 
time. From these plots, canopy pressure, strain and load pro- 
files will be prepared at times corresponding to those selected 
for photographic analysis. Pressure profiles will be input into 
the CANO program and the program pr dictions will be compared 
with the measured strain, shape and load data. 
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6.7.3 Flight Test Phase 

Preliminary Considerations Because of the cost and complexity 
of flight testing, the tests will be limited to low altitude 
>2500 ft) single parachute drops of one-third scale main para- 
chute models and a small number of full scale main parachute* 
drcps for verification purposes at high altitude (10,750 ft). 

This Is considered reasonable in view of the fact that the test 
Instrumentation and techniques involved are new and that 
interpretation of the data may be difficult, even for the 
simplest test conditions. Good profile shape data are considered 
to be of high value in correlating predictions of the CANO stress 
analysis program. Therefore, tes'.ing at low altitude to 
allow acquistlon of high quality photographs is considered a 
worthwhile compromise. 

The test objectives are summarized as follows: 

1 To obtain inflight measurements of 
riser force, pressure, strain and load 
distribution for one-third scale and 
full scale Apollo main parachutes, and 
to record photographically their shape 
and filling characteristics. 

2) To provide pressure distribution data 

for use in the CANO stress analysis program, 
and to provide strain, load and shape data 
for confirmation of program results. 


* Full scale main parachutes suitable for testing are presently 
available at Northrop Ventura. 
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3) To obtain data for substantiating the scaling 
laws developed in Reference 1, Section 2. 

Instrumentation The instrumentation will consist of: 

1) Pressure, strain and radial tape load transducers 
(whatever types appear most suitable after lab- 
oratory and wind tunnel tests), located in positions 
within the canopy identical to those of the winr 
tunnel models. 

2) Riser force, reefing and vent line load links. 

3) Ground-based motion picture cameras of focal length 
such that good resolution is obtained. 

4) Onboard cameras to record the parachute axial view 
during inflation. 

5) Airborne motion picture camera coverage for the full 
scale drops. 

All transducer outputs will be recorded by an oscillograph 
carried aboard the drop *-est vehicle, and dynamic pressure will 
be measured by an onboard pitot tube. All instrumentation wiring 
will be routed along radials and down suspension lines to the 
test vehicle. The data presentation format will follow that 
used in the wind tunnel tests. 

Test Conditions The scaling laws derived in Reference 1 for 
velocities and masses are, respectively: 

v^/v Q = (r^/r Q ) 5 (duplicates Froude number) 

M, /M q = (p^/p o ) (r^/r Q (duplicates added mass ratio) 

where the subscript 1 denotes the model and 
subscript o the full scale parachute 


217 


NVR 6432 



NORTHROP 

These expressions can be used to compute test conditions which 
make a model parachute simulate a full scale parachute. Table 
15 presents the required velocities and vehicle masses for the 
one-third scale parachutes for testing at 2500 ft altitude, 
simulating a main parachute opening at 10,750 ft after a high 
altitude abort. 


Table 15 . 

Test Conditions for 

Simulation 

of a 


Full Seal 

e Parachute by Use of 

a One- 



Third Scale 

Model 


Parachute 

Scale 

Altitude 

Velocity 

Vehicle Mass 



(ftl 

(ft/sec ) 

(lb) 

Apollo Main 

Fu]l 

10,750 

330 

6500 

Model 

1/3 

2,500 

190 

313 


Modeling the structural stiffness for dynamic simultude presents 
a serious problem. This is readily apparent if it is recalled 
that the strain is proportional to load per unit area (assuming 
linear strain curves for simplicity), while the load is propor- 
tional to the radius cubed. 


l£Zr 2 )- 

(F/r ; 


,2 \ 


~ 1L 


r p 
o o 


Thus the strains in the structural members of a geometrically 
scaled model (length, width, and thickness) would be only one - 
third that of a full scale model tested at the same air density. 
The model discussed in "ection 6.7.2 for wind tunnel testing (full 
scale thickness) would experience strains of one-ninth those 
of the full scale model, and the actual change in length (in 
inches) for a given member would be only 1/2/ the full scale 
value. It is obvious that the measurement of these small strains 
in fabric materials would present a problem. The strain distribu- 
tion, and therefore the shape, of the model would be slightly 
different from the full scale parachute because of the Improper 
stiffness modeling. 
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If, on the other hand, the construction details of a model 
with proper stiffness for dynamic simultude are considered, the 
difficulty of this approach is apparent. A one-third scale model 
of the Apollo main parachute with 68 gores would have suspension 
lines of 24 pound cord and radials of doubled 13 pound tape. 

These tapes would be equivalent to a strip of 1.1 ounce cloth 
(42 lb/in) with a width of 0.31 inches. This problem could be 
aleviated somewhat by using a smaller number of gores. For 
the sail material, however, it would be necessary to use some 
material other than nylon cloth to simulate the reduced stiffness 
while maintaining the proper porosity. 

For full scale parachute tests, it is considered that the most 
useful data can be obtained by testing at the actual flight 
altitude avid velocity of 10,750 feet and 330 ft/sec., respectively. 
The dynamic pressure at this condition is 90 psf and will require 
use of a pressure transducer with twice the range of the cur- 
rently available Melzlg-type (Schjeldahl) transducers. 

It should be noted that this limited test series Is designed 
to determine the suitability of the canopy instrumentation as 
well as to obtain useful data. It is anticipated that an ex- 
panded test series including drogue chutes and clustered para- 
chutes will be recommended as a next step. 


219 


NVR 6432 



SECTION 7.0 
SUMMARY 


The results of a one-year study of parachute structural analysis 
methods utilizing test data accumulated during the Apollo de- 
velopment and qualification test programs are presented., Study 
results include: 1) a literature review, 2) refinement and 

extension of the analysis methods, 3) corroboration of the 
methods by comparison of analytical and test results, 4) ap- 
plication of the Improved method to the Apollo parachutes, 

5) a study of dynamic loading effects in pilot parachute risers, 
and b) a study of the techniques of measuring loads, strains 
and differential pressures in parachutes. 

A method for computing the internal load distribution for a 
parachute with a given applied riser load and canopy differential 
pressure is presented. This analysis determines the canopy shape 
and Internal load distribution that satisfies equilibrium and 
boundary conditions for a given parachute. The validity of the 
analysis is demonstrated by comparing the analytically predicted 
shape with photographs of the Apollo parachutes. While man;, 
difficulties are encountered in the interpretation of film data, 
reasonably good correlation is obtained in comparisons of pre- 
dicted and observed shapes of the drogue, pilot and main para- 
chutes. (See Figures 21, 22, 32, 38, 39 » 4l and 42). 

Further verification of the method is obtained by comparing 
the predicted failure loads and modes for the drogue chute 
and three versions of the main parachutes with structural 
failures that occurred during the Apollo development testing. 

It is shown that the analysis correctly predicts the approxi- 
mate failure loads and location of the failure. (See Tables 1 
and 2) 
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The analysis Is Implemented in a computer program, CANO, which 
makes it possible to use a large digital computer f IBM 360/65 ' 
to perform the many iterative computations required for a solution. 
Use is made of this tool in performing several studies, as dis- 
cussed in the paragraphs below. 

Drogue Chute Structural Model Variation 

A study of the sensitivity of the internal load distri- 
bution to variations in radial tape stiffness was made 
to evaluate the assumptions that must be made concerning 
the effect of vertical members in the modeling of ribbon 
parachutes. Results show variations of 8 percent in 
horizontal member loading for the range of radial tape 
stiffness investigated. 

Drogue Chute Pressure Distribution 

Canopy profiles and internal load distributions are com- 
puted for four different pressure distributions for the 
full open drogue chute. For the two extremes investigated, 
i.e., peak pressure at the skirt and peak pressure at the 
vent, maximum horizontal ribbon loadir^ varies by 15 per- 
cent and maximum radial tape loadlrg varies by 10 percent. 

Drogue Chute Growth Study 

Results of a study of the variation of the shape and in- 
ternal loading in the drogue chute for different riser 
loading (with similar pressure distributions) show that 
the changes in diameter and internal loading are es- 
sentially linear. 
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Pilot Chute Internal Load Distribution 


Internal loads are computed for the pilot chute using 
a structural model that simulates the slots and the 
seams along the edges of the horizontal panels. Results 
show that the midsections of the center panels are the 
most critically loaded areas. 

Main Parachute Structural Models 

Three versions of the main parachute, representing major 
steps in the evolution of the present design, are com- 
pared. Results show that leading edge fullness in the 
lower sails causes high loadir^ In the trailing edges 
of these sails. Addition of reinforcing tapes to the 
trailing edges in critical areas lowers the load-to- 
strength ratio locally and also reduces the loading In 
the unreinforced areas. (See Figure 36 ) 

Main Parachute Failure Analysis 

Predicted failure loads and locations, based on the 
three structural models discussed above, are compared 
with structural failures that occurred during the Apollo 
development programs. Good agreement is shown between 
predicted and failure loads with a maximum spread of 
+ 13 percent. (See Table 2) 

Main Parachute Pressure Distribution and Film Analysis 

Since no information is available on the pressure distri- 
bution of reefed parachutes, the feasibility of inferring 
the pressure distribution from the shape of the canopy is 
investigated. Although the drop test films available do 
not give the profile views required for an accurate de- 
termination of pressure distribution, it is shown that a 
first approximation can be obtained from the data avail- 
able. This approximation should give better accuracy in 
the structural analysis than the assumption of uniform 
pressure used in the past. 
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Stress-Time Study oi’ the Main Farachute 

Previous analyses considered only three Instants lr. 
the opening process of the main parachute. It was 
assumed that maximum riser loading and maximum in- 
ternal loading occurred at the time that the parachute 
reached its maximum* inflation condition for the scage 
being investigated. In this study, drop test films 
are correlated with load traces to determine the riser 
load ,*.nd degree of inflation at eleven instants during 
the opening process. Pressure distributions are also 
estimated at these instants by the technique discussed 
in the preceding paragraph. With these inputs, program 
CANO is used to compute the load in each parachute ele- 
ment. From these results, stress -time histories for 
each sail and for a typical radial tape are constructed. 
(See Figures 45 through 48 and Table 4'. 

Optimum Weight Calculation for the Main parachute 


The weight of an idealized main parachute in which all 
members are Just strong enough at every point to with- 
stand design loading is given. If realistic Joint 
efficiencies are used In this calculation, the ratio of 
the Apollo main parachute weight to the optimum weight 
is 1.53. If 100 percent Joint efficiencies are used, 
this ratio is 2.17. Although these optimum weights 
can never be achieved in practice, they car. serve as 
a yardstick for Judging the efficiency of a parachute 
design and, together with proper statistical data, as 
a sound basis for predicting the weight of a new design. 
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Margin of Safety Calculations for the Main Parachute 

The computer program, CANO, that resulted from this 
study is applied in computing the margins of safety 
for the Apollo main parachute canopy to illustrate 
the reflneuients in structural analysis made possible 
by this study. The margins of safety given in 
Appendix B are the minimum values from 102 computer 
runs simulating the deployment process. The maximum 
applied riser loads for the one-drogue, two -main 
entry case are taken from Appendix C of Vol" 

Load-time histories are estimated from drc 
photographs. A comprehensive set of pressu 
tribution curves are applied at each openlrg *rcre- 
ment to cover the uncertainty in this area. 

A new internal loads analysis for ribbon paractu s is also 
included. In this analysis, called CANO 1, the structural model 
includes vertical ribbons. Internal load solutions are given 
for the drogue chute In both the reefed and disreefed condition. 
The results show that a significant portion of the meridional 
load Is carried by the vertical ribbons, thereby giving lower 
radial tape loading. No tust data to corroborate these results 
are available, so the analysis must be used with caution. 

A study of data obtained in dynamic loading tests of pilot chute 
risers is given. An approximate linear-elastic analysis produces 
explanations of many aspects of the test data.' It is shown that 
longitudinal wave effects must be included in the structural 
analysis of pilot chute risers since multiple failures and 
failures in free lengths may occur at loads that would not cause 
failure if applied statically. It is also shown that hystersis 
effects must be considered in the load-strain relationships used 
in the riser dynamic analysis. Improved results are obtained by 
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use of nonlinear static load-strain curves and a more accurate 
representation of the boundary conditions in a computerized finite 
difference method approach. Further refinement of the mathemati- 
cal model is shown to be necessary. Analysis of the pilot chste 
riser test data shows that the dynamic load-strain curves followed 
by the riser nylon webbing are approximately linear, but that 
small variations from linearity produce effects which cannot be 
neglected. It is seen that the longitudinal wave speed is 
sensitive to some as yet undiscovered variable in the tests 

With the facility provided by the CAHO analysis for detailed 
structural modeling and for investigating many loading conditions, 
the accuracy of the reported structural analysis is limited 
mainly by the accuracy of input loading data. Although accurate 
methods for predicting riser loads have been developed (see 
Volume I), little information is available on the distribution of 
aerodynamic forces on the canopy. It Is shown In the analysis 
that this distribution has a first order effect on the loading 
of individual parachute members. Results of a study of the 
technique 0 of measuring canopy differential pressure distribution 
and Internal loading are given, together with a plan for a three- 
phase test program to acquire these data for the Apollo parachutes. 
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SECTION 8 
CONCLUSIONS 


The conclusions based on the work described in this report are 
as follows: 

1) The structural analysis method Implemented in 
computer program CANO gives a satisfactory esti- 
mate of the Internal load distribution in a para- 
chute under known riser and aerodynamic loading. 

This conclusion is based on comparisons of pre- 
dicted shapes with drop test photographs and on 
comparisons of failure loads and modes with test 
results. 

2) When good photographic data and riser load traces are 
available, it is possible to infer the canopy 
differential pressure distribution from the canopy 
shape. Typical drop test photographs do not pro- 
vide the desired accuracy, however. 

3) The analysis techniques used in the stress-time 
study of the main parachute can be used to pre- 
dict the state of stress in a parachute throughout 
the opening process provided that accurate shape-time 
and riser force history data are available. 

4) An analysis of the Apollo drogue chute (a conical 
ribbon parachute) indicates that a significant 
portion of the meridional load is carried by the 
vertical members. Corroborating test data, which 
could be obtained by measuring vertical ribbon and 
radial tape loads, are needed to verify the analysis 
presented . 
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5) Most of the unexpected results observed In 
the riser dynamics tests can be explained by 
a graphical linear-elastic analysis. A more 
accurate, nonlinear-elastic analysis can be made 
by using finite difference methods to solve the 
dynamical equations of parachute risers mder 
Impact loading. 

6 ) Unusual failures such as simultaneous multiple 
breaks and failures in free lengths of pilot 
chute risers can be caused by tensile strain 
waves at loads that would not cause failure if 
applied statically. The wave speeds observed 
in the tests were sensitive to some as yet un- 
determined variables. 

7) The most important unknown factor in parachute 
structural analysis is canopy differential pres- 
sure distribution. A method for predicting pres- 
sure distribution is needed for pretest structural 
analysis. 
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SECTION 9.0 
RECOMMENDATIONS 

Based on the results of the one-year studj of structural 
analysis techniques for Apollo spacecraft parachutes, it i 3 
recommended that: 

1) The stress-time analysis method presented in 
this report be used in future parachute design 
and development programs. 

2) The internal load analysis, CANO, be used to make 
a systematic study of the effects of the major 
design parameters on parachute weights. 

3) The concept of a theoretically optimum parachute 
be used for evaluating the structural efficiency 
of existing parachute designs, and that it be used 
as the basis for predicting the weight of future 
designs. 

4) The study of the variation of canopy shape with 
pressure distribution be continued to provide a 
method for estimating pressure distribution from 
drop test photographic data. 

5) A test program be undertaken to measure pressure, 
strain and load distribution in parachutes during 
the opening process. 

6) Corroboration be sought, by test measurements, for 
the internal load analysis for ribbon parachutes 
developed in this report. 

7) A method be developed for predicting differential 
pressure distribution on parachute canopies, pre- 
ferably an analytically based method with empirical 
support; this is required to remove a barrier to 
advancement of parachute stress analysis technology. 
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8) The usefulness of photographic data be improved 

In future flight tests by: a) use of airborne, 

high resolution cameras to obtain good profile views 
at the critical stages of each test, b) distinctive 
marking of the canopy members, and c) calibration 
of onboard cameras before each flight. 

9) The variable wave velocity and other phenomena be 
investigated by conducting tests similar to those 
described herein, but with greater control over the 
manner in which the Impact load is applied to the 
end of the specimen. The specimens should be de- 
signed in such a way that discontinuities do not 
Interfere with precise observation of the wave 
action. One series of tests should be optimized 
with respect to determining material properties 
(dynamic load-strain curve) and a second series of 
tests optimized with respect to determining the 
behavior (transmission and reflection coefficients) 
of joints, links and other riser components. 

10) The AFTON computer program be applied to a complete 
riser model using material property and component 
performance data as determined in (9). 

11) Experimental Impact studies be conducted using 
complete riser system specimens in order that the 
effect of the true boundary conditions be as- 
certained. 

12) That the data obtained in (9), (10) and (11) be 
used to develop a qualified method for dynamic 
analysis of pilot parachute risers. 
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APPENDIX A 


SUMMARY OF APOLLO AERIAL DROP TESTS 


The following table contains a summary of the Apollo aerial drop 
tests that have been conducted to date. This Information can be 
used as a quick reference to find the test objectives, type of 
vehicle, parachute versions, and recorded loads for any particular 
drop test. The listed references indicate where this data was 
obtained and where additional data may be found. 

Three types of reports are used to record the drop test Information: 
Paradynamics Technical Memoranda (PTM), Flight Test Reports (FTR), 
and Northrop Ventura Reports (N7R). 

The Paradynamics Technical Memorandums and the Northrop Ventura 
Reports are formalized reports which were used to compile the in- 
formation resulting from the Apollo drop test program. The Para- 
dynamics Systems Section was responsible for preparing the PTM's, 
and the Field Test Branch was in charge of the preparation of the 
NVR's. Flight Test Reports are data packages which contain test 
objectives, test data, and post-test evaluation of aerial drop 
tests. These reports were, prepared by the Analytical Engineering 
Branch. 

Reference is sometimes made to tests as being part of Block I, 

Block II or Block II (H). The tests associated with these blocks 
are as follows: 

Block I Test Series 1 thru 62 and 71 

Block II Test Series 70 and 73 

Block II (H) Test Series 80 thru 99 
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Table A-l (Continued) 
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Table A-l (Continued) 
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APPENDIX B 


EXAMPLE MARGIN OF SAFETY CALCULATIONS FOR THE APOLLO 

MAIN PARACHUTE 

Margin oY safety calculations for one design condition for the 
Apollo main parachute, and a comparison with the values given 
in the Block II (H) stress analysis'^ are presented in this 
appendix to illustrate the refinements in parachute structural 
analysis made possible by this study. 

Design loads are based on the maximum loads for the one-drogue, 
two-main (1D/2M) entry condition computed in Appendix C of 
Volume I. * A load-time history corresponding to these maximum 
stage loads is constructed analytically with the aid of drop 
test films and load traces from representative flight conditions 
and configuration. Canopy differential pressure distributions 
are estimated from shape data observed in drop test films and 
from the wind tunnel data of Reference 20. The uncertainlty 
in this parameter is satisfactorily resolved by applying a set 
of pressure distribution curves that covers the range that can 
be reasonably expected to occur. Loadings for all elements of 
the canopy structural model (see Section 4.3.1) are computed by 
means of program CANO. Safety factors, temperature degradation 
factors, and joint efficiencies are taken from the Block II (H) 
stress analysis. Reference 19> PP 281-292. 

RISER LOADING 

The design ultimate loads for the one-drogue, two-main entry 

case are obtained by multiplying the limit loads computed in 

this study (Appendix C of Volume I) by a safety factor of 

1.35. The table below gives the new loads used in tnis analysis, 

and also the corresponding values, for comparison, from Reference 19. 
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Design Ultimate Loads - 1D/2M Entry Condition 


First stage 
Second Stage 
Open 


^Cev^_Loac^ 

(1.35 x Limit Load 
from Volume I) 

25,200 lb 
24,800 lb 
25,200 lb 


Old Load 

(From Reference 10) 

27,350 lb 
31,6oo lb 
29,700 lb 


For this analysis. Internal loads are computed at critical time 
Instants throughout the opening process. These events, called 
inflation states in this analysis, correspond to the inflation of 
the Individual sails. For example. Inflation State 1 occurs just 
as Sail 1 inflates and ceases to flutter randomly. By identifying 
this event in drop test films and correlating the t:\ming with 
load traces, the fraction of the maximum stage load present at 
each inflation state can be established. Design loads are ob- 
tained by multiplying the stage design load by this factor. 

Table B-l summarizes the design loads used in this analysis. 


Inflation 

Table B-l 
Fraction of 

Design 

State 

Max. Stage 

Ultimate 


Load 

Load, lb 

1 

.741 

18,670 

2 

. 806 

20,310 

3 

.882 

22,200 

4 

1.00 

25,200 

5 

.819 

20,310 

6 

1.00 

24,800 

7 

.819 

24,800* 

8 

• 591 

24,800* 

9 

. 645 

16,250 

10 

.795 

20,030 

11 

1.000 

25,200 

12 

.952 

23,990 


* Maximum load is applied throughout latter portion of this stage 
to cover uncertain! ty »n inflation-time history. 
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PRESSURE DISTRIBUTION 


Program CANO proportions local pressure acting on each element 

in accordance with any sDecifled distribution curve and adjusts 

overall magnitude to balance the applied riser load. Variations 

in the shape of the pressure curve have a first order effect on 

the stresses in individual elements. Because of the many practical 

problems involved in direct measurement of differential pressures 

in irflatlng parachutes, no measurements were made during the 

Apollo test programs. Some knowledge of this parameter has been 

gained indirectly from drop test photographic coverage by comparing 

observed canopy shapes with shapes computed by the CANO computer 

program, using assumed distribution curves (see Section 4.3.3^* 

20 21 

Fressure measurements made by Melzig * in wind tunnel and 
flight tests, while not directly applicable to the Apollo para- 
chutes, provide general knowledge of the variations in pressure 
distribution during Inflation. However, the available data is 
not considered accurate enough for prediction of precise pressure 
distributions for use in the stress analysis. To provide a con- 
servative analysis, a comprenensive set of distribution curves 
that covers the worst cases that can be reasonably expected to 
occur is applied for each inflation state. xhe family of curves 

shown in Figure B-l resulted from shape analysis of Apollo drop 

20 21 

test films, from Melzig' s test results * , and from a computer 

study of the effect of pressure distribution curve shapes on in- 
ternal loading. As shown in Figure Bl, a pressure peak is applied 
at various locations in the inflated portion of the canopy from 
the vent to the leading edge. For each inflation state, load 
calculations for the entire canopy are made for a number of pres- 
sure curves, selected to give peak pressures on each sail in turn. 
Table B-2 gives a summary of the results of these calculations in 
the form of load -strength ratios (applied load t rated strength) 
for the critical element of each sail for the loading conditions 
considered. 
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Load 

lb 

Pressure 
[ Curve No. 

First Stage 

18,670 

9 


Table B-2 Sail Load-Strength Ratios 


Sail Number 



20,310 



.36 


* .37 .07 


.34 .18 


, .36 .25 .03 


22,220 

16 

.38 .33 .04 

22,220 

3 

.35;. 34 .10 

25,200 

9 

•40 1 .35 .06 .02 1 : : 

j - ■ ’ .. i i 

l 

^ 17 

.39 .38 .08 .02 


5 

.37,. 37 .19.03 

\ 

12 

•35. 35 .261.06 

.... - ..... - I ■ T~ - 1111 - - - m - + ,, (ll , , m , ■ 

25,200 

20 

.32;. 33 .28. 10 


Second Stage 


20,310 9 

1 8 


' | — 

20,310 ! 10 


.31: . 30 .10, .02 0 

T731 -31 .13^.02 0 _ __ 

.29 .30 .20—03. 0 
~ .27 .28 .23 .07 0 ' 

.25 .26 .23 .09 .04 


I -35 .34 .20 .04 .02 .02 


I .34! .34 .22 .04 .02 .02 


.33 .34 .26 .06 .02 .02 
.29;. 31 .28 .12 .03 .02 


.27 * .29 .28 .19 .09 .02 


.27 .28 .28 .2 0 .11 .09 
.33]. 33 .2^ . 07 .03 j 03 -of 
^32 .33 .03 .03 _.0l 

.30^.32 .28 .12^.04 . 03 .01 


.28:. 30 .29 .18 .05 .04 .02 


.251.28 .28 .23 .12 . 10 .0 
.23 >.26: .27 .24^.16 .15 .08 
31 1 .327 .26^ . 127.04 .03 .0 

0 


6 

.27 

j -29 j 

.29; .22 
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Table B-2 (Concluded) 


Load Pressure 
lb Curve No . 


24,600 



24,800 


Third Stage 


16,250 9^ 

8 _ 

16 

* 14 


16,250 


20,030 


20,030 


25,200 

! 


25,200 


23,990 




23,990 


.19 


.17 .16 .19 


.14 .16 .17 


.13 .14 .16 


.12 .14 .15 


.12 .13 


.16 


.14 


.14 .15 .19 


.13 .14 .16 


.06 .12 .15 


.23 .25 


.20 .23 .25 


.17 .20 .23 


.20 


.13!. 16 .19 


.13 .15 .18 


.12 .14 .17 


.20 .23 .24 


.12 .19 .22 


.15 .17 


.12 


.04 .10 .13 


.03 


.04 


.16 .06 .05 


.16 .13 .13 


.16 .16 .08 


.15 .18 .21 


.15 .19 .22 


.19 


.22 .25 


.18 


.17 .24' .2 8 



.291 

.33 

Os 

cv 

• 

.34 

.29 

.34 

.29 

.33 

to 

CM 

• 

.33 

.19 

.21 

.27 

A 1 

.29 

.34 

.26 

.33 

.23 

.29 


.02 .01 


.01 0 


.02 


.04 


.16 .07 


.19 .10 


.12 .04 


.19 .06 


.24 .12 


.27 .18 


.08 .04 


.13 .05 


.22 .09 


.29 .15 


.33 .23 


.26 .29 


.35 .30 


.35 .30 


.06 


.15 


.21 


.35 .32 


.33 .35 


- 4 - 



.04 


.02 .01 


.03 .01 


.07 .02 


.10 .05 


.02 . 01 | . 01 ; 


.03 


.04 


.05 


.23 .13; .06 


.03 


.05 


.08 .03 


.21 .09 .04 .01 


.31 .28 
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MARGINS OP SAFETY WITH NEW LOADS 

The critical load-strength ratios from the CANO solutions are 
used, together with Joint efficiencies and design factors from 
the Block Il(H) stress analysis^, to compute the margins of 
safety for the major canopy components. 

j-g _ Joint Eff x Thermal Factor x Abrasion Factor _ . 
Load-Strength Ratio x Unsym. Load Factor 

Corresponding margins from Reference 19 are also shown for 

comparison. The difference in the compared margins is due to 

the revised load input values and the refined stress analysis. 

Sails 

The critical load-strength ratios shown in Table B-2 are used, 
together with the design factors listed below, to calculate the 
margins shown in Table B-3. 

Design factors: 


Joint Efficiency 

0.55 

Thermal Factors 


Abrasion Factor 

1.C0 

First Stage 

0.87 

Unsym. Loading 

1.05 

Second Stage 

0.92 



Open 

0.94 

Table B-3* Sail Margins 

of Safety with New Loads 

Sail Number 


Margin of Safety 


1st Stage 

2nd Stage 

Open 

1 

.15 

.48 

1.5 

2 

.22 (.05)* .40 

1.2 

3 

.63 

.67 (.17)* 

1.0 

4 

high 

.90 

1.05 

5 


.82 

.69 

6 


.89 

.46 

7 


1.6 

.40 

8 


high 

.43 

9 



.58 

10 



CVJ 

• 

11 



1.2 

12 



high 


* Values shown in parenthesis ( ) are margins givon in Reference 
19 for the 1D/2M entry condition. 
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Suspension Lines 


Maximum suspension line loads from the calculations dicussed 
above are: 

First stage 371 lb 

Second stage 367 lb 

Third stage 379 lb 

Rated strength and design factors from Reference 19 are: 


Rated strength 

650 lb 


Joint efficiency 
Thermal factors 

0.86 


First stage 

0.87 


Second stage 

0.92 


| Open 

0.94 


Abrasion factor 

0.96 


Unsymmetrical loading 1.05 


Margins of Safety with New 

Loads 

Present Analysis 

From Reference 19 

First stage 

0.20 


Second stage 

0.28 

0.0 

Open 

0.27 


Skirt and Vent Bands 
Maximum loads computed for 

the skirt and vent 

bands for the full 

open, split gore condition 

are: 


Skirt band 

1240 lb 


Vent band 

3590 lb 


Rated strength and design factors from Reference 19 are: 

Rated strength Skirt band 2400 1. 



Vent band 5200 . 
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Joint efficiency o,95 

Abrasion o . 9>' : * 

Thermal 0.96 

Margins of Safety with New Loads 

Skirt band = 0.69 

Vent band = 0.27 

Corresponding values are not shown in Reference 1? since che 
entry condition was not crlcital for these members. 

Radial Tapes 

Maximum radial tape loads from the previous calculations are: 


First stage 

381 lb 


Second stage 

415 lb 


Open 

379 lb 


Rated strength and design 

factors from Reference 

19 are: 

Rated strength 

2 x 390 = 700 lb 


Joint efficiency 

0.89 


Abrasion factor 

0.96 


Unsymmetrical load 

1.05 


Thermal 



First stage 

O.87 


Second stage 

0.92 


Open 

0.96 


Margins of Safety with New 

Loads 


present Analysis From 

Reference 19 

First stage 

MS * .30 


Second 3tage 

MS = .26 

.09 

Open 

MS = .44 
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CONCLUDING REMARKS 


This analysis could be extended to other design cases where 
applicable empirical data (load-inflation history) is available. 
It Is of Interest to note here that the calculations for each 
of the 68 solutions reported in Table B-2 would have required 
approximately one man-week of effort with the analysis tools 
available at the beginning of this study. With program CANO , 
additional solutions (once the structural model is set up) are 
obtained by a single punched card input. Execution time on a 
large digital computer (IBM 360/65) averages less than two 
minutes per solution. 

The example case (one-drogue, two-main, entry) is presented to 
illustrate the Improved stress analysis made possible by this 
study. The foi*egolng MS calculations are therefore intended to 
provide an illustration of the techniques and not to establish 
official changes in the Apollo ELS margins of safety of record. 
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